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Abstract 

A  Code  Division  Multiple  Access  (CDMA)  spread  spectrum  packet  radio  network 
model  is  presented.  The  topology  is  a  single-hop  fully  connected  network  with  identi¬ 
cal  users.  The  network  model  allows  the  performance  of  the  radio  links  to  be  specified  in 
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into  account  the  restriction  that  half-duplex  radios  cannot  transmit  and  receive  simultane¬ 
ously.  Single  hop  throughput  and  probability  of  successful  first  transmission  are  derived. 
The  effects  on  throughput  of  the  received  signal  strength  E^/Nq,  the  number  of  chips 
per  bit  N  and  the  number  of  users  M  are  shown.  A  channel  load  sense  access  protocol 
is  introduced  in  which  radios  are  blocked  from  transmitting  when  the  channel  is  heavily 
loaded.  The  increase  in  throughput  due  to  this  protocol  is  shown  for  zero  and  for  non-zero 
propagation  delay. 
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case  for  half-duplex  radios. 

Earlier  CDMA  network  models  assume  a  network  access  protocol  which  is  analogous 
to  ALOHA;  namely,  if  a  packet  arrives  at  an  idle  radio  it  is  immediately  transmitted.  For 
ALOHA  systems,  network  performance  can  usually  be  improved  by  introducing  carrier 
sensing,  where  the  users  listen  before  transmitting  and  block  transmission  if  another  carrier 
is  sensed.  In  this  report,  we  investigate  the  improvement  that  is  realized  by  using  a  channel 
load  sense  access  protocol  in  a  CDMA  network.  In  this  protocol,  radios  are  allowed  to  begin 
transmitting  a  packet  only  if  the  total  number  of  transmissions  on  the  channel  is  less  than 
a  threshold  K. 

We  analyze  specific  fully  connected  configurations  of  a  more  general  multi-hop  model 
developed  by  Boorstyn  and  Kershenbaum  [BOOR80]  and  also  by  Brazio  and  Tobagi 
[BRAZ84],  For  these  configurations,  we  find  the  single-hop  throughput  and  the  proba¬ 
bility  of  correct  packet  reception.  The  models  are  for  asynchronous  networks  with  variable 
length  packets  that  are  generated  by  Poisson  processes.  In  Section  2,  we  present  the  sys¬ 
tem  assumptions  underlying  the  models  analyzed.  In  Section  3,  we  present  a  simplified 
network  model,  which  assumes  a  simplified  channel  model.  We  analyze  the  throughput 
and  the  probability  of  a  packet  being  successful,  and  present  some  performance  results. 
In  Section  4,  we  present  a  more  accurate  channel  model  which  can  include  forward  error 
correction  (FEC)  coding.  We  refine  the  network  model  to  accomodate  this  channel  model. 
In  Section  5,  we  modify  the  network  model  in  order  to  investigate  the  performance  of  a 
network  of  half-duplex  radios.  Then,  in  Section  6,  we  add  channel  load  sensing  to  the 
models  developed  earlier,  and  investigate  the  resulting  improvement  in  performance.  By 
using  an  approximate  model  of  propagation  delay,  we  find  the  degradation  in  throughput 
that  this  delay  causes  in  networks  using  channel  load  sensing. 

2.  System  Assumptions 

The  general  model  applies  to  various  networks  with  the  following  common  features. 
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1.  Introduction 


For  many  years,  various  forms  of  spread  spectrum  signalling  have  been  used  in  secure 
communications  systems.  Though  primarily  implemented  in  military  systems  because  of 
their  anti-jam  and  anti-intercept  properties,  spread  spectrum  techniques  offer  other  ca¬ 
pabilities  such  as  ranging  and  code  division  multiple  access  (CDMA).  The  performance 
of  CDMA  systems  which  use  forward  error  correction  (FEC)  coding  has  been  analyzed 
by  various  authors.  Raychaudhuri  [RAYC81]  analyzed  throughput,  delay,  number  of  re¬ 
transmissions  and  stability  of  a  slotted  single  hop  network.  Musser  and  Daigle  [MUSS82] 
derived  the  throughput  for  a  single  hop  unslotted  network  with  fixed  length  packets.  Purs- 
ley  [PURS83]  studied  asynchronous  frequency  hopping  systems  with  fixed  length  packets, 
to  find  the  throughput  and  the  probability  of  a  packet  being  correctly  received.  Taiple 
[TAIP84]  investigated  the  channel  and  network  performance  of  a  Direct  Sequence  CDMA 
(DS-CDMA)  packet  radio  system. 

Previous  models  and  analyses  of  the  performance  of  CDMA  systems  fall  into  one  of 
two  groups.  The  first  group  focuses  on  the  radio  links,  deriving  bit  error  rates  (BER)  for  a 
given  code  and  number  of  interfering  transmissions.  The  other  group  focuses  on  network 
performance,  deriving  throughput  and  delay.  Typically,  in  the  analyses  of  network  perfor¬ 
mance,  a  very  simple  channel  model  is  used  to  approximate  the  radio  link  performance. 
Also,  the  models  only  account  for  the  transmitters,  assuming  that  a  receiver  is  always 
available. 

In  this  report,  we  present  a  network  model  which  is  flexible  enough  to  incorporate 
a  variety  of  channel  models  for  the  radio  links.  We  than  adopt  a  channel  model  of  a 
Direct-Sequence  Binary  Phase  Shift  Keying  (DS-BPSK)  radio  link  presented  by  Taiple 
[TAIP84],  and  find  the  network  throughput.  A  modification  of  the  network  model  allows 
us  to  account  for  the  receivers  as  well  as  the  transmitters.  We  investigate  the  throughput 
of  a  network  in  which  radios  cannot  transmit  and  receive  simultaneously,  as  would  be  the 
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The  network  is  fully  connected  with  only  single  hop  traffic.  All  users  are  identical.  In 
Section  5,  we  use  an  M  user  finite  population  model,  and  consider  each  user  to  have  a 
directed  link  to  each  of  the  M  —  1  other  users.  We  assume  equal  traffic  requirements  on 
all  M(M  —  1)  of  these  links. 

The  network  has  the  following  characteristics  due  to  the  use  of  spread  spectrum  CDMA. 
Radios  are  assumed  to  use  receiver  directed  codes.  The  destination  receiver  has  a  probabil¬ 
ity  a{  of  synchronizing  to  the  packet  header  (capturing  the  packet),  where  at  is  a  function 
of  the  number  of  radios  transmitting  at  the  start  of  the  packet  transmission.  Once  cap¬ 
ture  occurs,  we  assume  that  the  receiver  never  loses  synchronization  until  the  packet  is 
completed. 

Two  channel  models  are  considered.  The  first  is  a  simple  model  in  which  no  trans¬ 
mission  errors  occur  for  L  or  fewer  simultaneous  transmissions,  and  transmission  errors 
occur  with  probability  1  for  greater  than  L+  1  simultaneous  transmissions.  This  is  referred 
to  as  the  step  function  channel  model,  as  the  probability  of  bit  error  vs.  the  number  of 
interfering  transmissions  is  a  step  function. 

The  second  channel  model  is  a  more  accurate  model  of  a  DS-BPSK  CDMA  system. 
Channel  performance  results  are  found  for  a  system  incorporating  FEC  coding  as  well 
as  for  an  uncoded  system.  The  type  of  coding  modeled  is  convolutional  coding  with 
Viterbi  decoding.  In  our  analysis,  we  assume  an  upper  bound  L,  such  that  more  than 
L  simultaneous  transmissions  will  cause  all  packets  to  be  unsuccessful  with  probability 
1.  Of  course,  there  is  a  small  probability  that  a  packet  will  be  successful  even  if  the 
number  of  interfering  transmissions  is  so  great  that  the  probability  of  bit  error  is  nearly 
1/2.  Nevertheless,  the  time  averaged  contribution  to  throughput  by  such  packets  is  very 
small.  We  choose  the  cutoff  point  L  such  that  the  error  due  to  ignoring  larger  numbers  of 
interfering  transmissions  is  negligible. 

In  both  channel  models,  we  assum®  that  a  single  error  in  the  corrected  information 
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bit  stream  will  cause  the  packet  to  be  unsuccessful.  In  order  to  identify  packets  with  bit 
errors,  parity  or  cyclic  check  bits  would  be  appended,  and  those  packets  found  to  be  in 
error  would  be  discarded.  We  will  ignore  the  overhead  this  parity  check  would  require. 

We  ignore  the  effect  of  acknowledgments,  assuming  a  perfect  and  instantaneous  ac¬ 
knowledgement  channel  is  available. 

A  packet  which  is  scheduled  for  transmission  is  blocked  (i.e.,  rescheduled)  if  the  radio 
is  currently  transmitting  or  receiving.  In  the  channel  load  sense  models,  the  transmission 
is  also  blocked  if  the  number  of  radios  currently  transmitting  is  greater  than  or  equal  to  a 
threshold  K.  Packets  that  are  not  blocked  are  immediately  transmitted  at  the  scheduling 
point 

For  channel  load  sense  models,  it  is  assumed  that  each  radio  has  perfect  knowledge  of 
the  number  of  radios  currently  transmitting.  At  a  scheduling  point,  the  decision  whether 
or  not  to  transmit  is  based  on  this  knowledge.  In  a  real  system,  the  users  can  only  estimate 
the  channel  loading.  For  example,  if  headers  are  transmitted  using  a  general  code  rather 
than  a  receiver  directed  i-ode,  the  user  can  count  the  Leaders  received  in  the  recent  past 
to  estimate  the  number  of  other  users  transmitting.  Alternatively,  the  received  signal  plus 
noise  power  can  be  integrated  to  give  an  estimate  of  the  total  energy  in  the  pseudo-noise 
signals  received,  which  will  indicate  the  channel  loading.  Clearly,  these  implementations 
will  perform  worse  than  the  idealized  model. 

3.  Network  Model 

In  the  basic  model,  we  assume  that  receivers  are  dedicated  to  transmitters,  so  that 
transmissions  are  never  blocked  due  to  the  source  being  busy  receiving,  and  transmis¬ 
sions  are  never  unsuccessful  due  to  the  destination  being  busy.  For  finite  transmitters, 
each  user  generates  scheduling  points  from  a  Poisson  process  at  the  rate  of  offered  traffic 
A.  This  Poisson  traffic  includes  retransmissions  ?s  v/eH  as  newly  generated  pa.kcts.  We 
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use  the  heavy  traffic  assumption  that  there  is  a  packet  available  for  transmission  at  ev¬ 
ery  scheduling  point.  If  the  user  is  not  already  transmitting  and  the  protocol  does  not 
inhibit  transmission  (in  the  channel  load  sense  protocol),  the  transmission  of  the  packet 
will  begin  at  the  scheduling  point.  For  an  infinite  population,  the  scheduling  points  for 
all  transmissions  are  generated  from  a  Poisson  process  with  aggregate  rate  A.  Again,  a 
packet  will  begin  to  be  transmitted  at  the  scheduling  point  if  the  protocol  does  not  inhibit 
transmission.  The  transmission  time  of  the  packets  is  exponentially  distributed  with  mean 
l//x.  We  use  g  to  denote  the  normalized  rate  of  offered  traffic,  g  =  A//z. 

In  this  simplified  model,  we  assume  a  channel  behavior  such  that  no  errors  occur 
for  L  or  fewer  radios  transmitting,  and  errors  occur  with  probability  1  for  more  than  L 
simultaneous  transmissions.  A  more  realistic  model  of  the  channel  behavior  is  presented 
in  Section  4. 


Let  X{t)  be  the  number  of  radios  transmitting  at  time  r  (We  use  r  to  denote  time 
throughout  the  report  to  avoid  confusion  when  we  use  t  as  the  number  of  active  transmitters 
in  Sections  5  and  6).  We  find  that  X(r)  is  a  continuous  time  Markov  chain.  The  state  space 
S  of  X(r)  is  5  —  {0,1,2,....}  for  an  infinite  population,  and  5  =  {0, 1,  2, ...,  M }  for  finite 
transmitters.  In  both  cases,  5  has  a  strictly  positive  stationary  distribution  {jry,  6  5}. 
We  notice  that  X(r)  for  the  infinite  population  model  has  the  same  state  transition  rates 
as  the  model  of  the  M/M/ oo  queue,  shown  in  Fig.  3.1,  and  for  the  finite  transmitters 
case,  the  same  as  the  M/M/oo/ /M  queue,  shown  in  Fig.  3.2  [KLET75].  For  an  infinite 
population, 


' 


j  =  0, 1,2, ... 


For  finite  transmitters, 


(3.1) 


_  _ (V __  (M\ 

1  (1+A  M»\j) 


j  =  0, 1, ...,  M 


(3.2) 
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Figure  3.1.  State-transition-rate  diagram  for  the  infinite  population  model. 


Figure  3.2.  State-transition-rate  diagram  for  the  finite  transmitters  model. 
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3.A  Probability  of  Saeeess 

We  derive  the  probability  of  success  for  the  infinite  population  model  by  tracing  the 
evolution  of  the  network  model  during  the  transmission  of  a  tagged  packet.  If  the  transmis¬ 
sion  is  completed  before  the  state  L+ 1  is  reached,  the  packet  will  be  successful.  Conversely, 
if  the  state  L  +  1  is  visited  first,  the  packet  will  be  unsuccessful. 

We  define  the  following  parameters: 

User  t  is  the  source  of  the  tagged  packet 

to  is  the  time  at  which  the  transmission  of  the  tagged  packet  begins 
Ps  is  the  probability  that  a  transmitted  packet  is  successful 

Ps ij  is  the  conditional  probability  that  a  transmitted  packet  is  successful  given  that 
the  network  was  in  state  j  at  time  Tq 

T}  is  the  steady  state  probability  of  the  network  being  in  state  j  at  time  r0~  given  that 
a  packet  transmission  begins  at  time  To 

Pin>lF\j)  is  the  probability  that  a  tagged  user  t  is  not  transmitting  given  that  the 
network  state  is  j 

To  find  Ps,  we  first  condition  on 

oo  L—  l 

ps  =  E  *)ps\j  =  E  *jps\j  (3-3) 

J=0  J= 0 

The  tagged  packet  began  transmission  at  time  to,  which  implies  that  user  t  was  idle  at 
time  Tq~  .  Thus, 

Kj  =  Pr(  j  transmitters  |  user  t  was  idle) 

For  an  infinite  population,  due  to  the  Poisson  generation  process,  iry  =  Try.  For  finite 
transmitters,  we  find  from  Bayes’  Theorem  as 

^  _  PxIDLfAi}* j 

1  r£n‘  Pi,DL'U)*i 
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Because  all  users  are  identical, 


M-j 

M 


for  t  =  1, 2, .  .,  M 


Thus,  we  have 


(3.5) 


L- 1 


Y  Kjps\j 


infinite  population 


E"o‘(Af  -  i)*jPs,, 

E"o‘(W-J>y 


finite  transmitters 


(3.6) 


If  we  count  the  number  of  successful  packets  transmitted  in  a  long  interval  of  time, 
and  divide  by  the  length  of  this  interval,  we  will  get  a  variable  whose  expected  value  is  7, 
the  number  of  successful  packets  over  time.  We  can  find  7  from  Psy  as 


E  Wsu 

J=0 

min(L,A/)  —  1 
J=d 


infinite  population 


finite  transmitters 


(3.7) 


For  the  simple  case  of  an  infinite  transmitter  population  with  L  =  1,  which  is  a  standard 
ALOHA  channel,  ttq  =  exp(  —  A//x)  and  Ps\0  is  the  probability  that  no  arrivals  occur  before 
the  packet  is  completed.  Both  the  time  to  completion  of  the  message  and  the  time  until 
the  next  arrival  are  exponentially  distributed,  so  we  can  easily  find 


(3.8) 


Ps  = 


e'9 

1  +  g 


(3.9) 
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and 


1  _  ge  g 

fi  1  +  g 


(3.10) 


This  case  was  previously  analyzed  by  Ferguson  [FERG75]  and  by  Bellini  and  Borgonovo 
[BELL80]. 


For  the  more  difficult  case  of  L  >  1,  we  can  find  PSy  recursively  as  follows.  A  packet 
finding  the  system  in  state  j  will  cause  a  transition  to  state  j  4-1.  From  state  j  +  1,  we 
have  the  following  probabilities  for  the  next  transition. 


PAj)  =  Prfnext  transition  is  due  to  an  arrival)  =  r - - 

K  (j  +  l)/x  +  A 

PB{j)  —  Pr(next  transition  is  due  to  another  packet  completing) 

__  _ W _ 

(j  +  l)H  +  \ 

Pc{j)  =  Pr(next  transition  is  due  to  the  tagged  packet  completing) 

_  _ t _ 

(;  +  i)a*  +  ^ 

Because  of  the  memoryless  property  of  the  exponential  packet  length  distribution,  at  every 
state  transition,  given  that  no  errors  have  occurred  and  that  the  state  entered  is  state  j, 
the  probability  that  the  packet  will  be  successful  is  the  same  as  Ps\)-\-  So  we  have  the 
equations 


(3.11) 

(3.12) 

(3.13) 


Ps\o  —  Pa{q)P s|i  +  ^c(O)  j  —  0 

=  PaU)Ps\j+i  +  PtU)Ps\j- 1  +  Pc(i)  1  <  j  <  L  —  1 
Ps\j  =  0  j>L 


0 


or 


0  =  -(/i  4  A)P5|0  4-  A P 5jj  4-  /i 

0  =  “((;  +  l)/i  +  A)P5|j  +  APS|J+i  +jnPS\j-i  +fi  \<j<L-\  (3-15) 

0  =  ~(L/l  +  A)P 5|£_x  +  (L  —  l)pPs\L-2  +  P 

We  define  •  to  be  the  vector  such  that  [•]_,  =  P5jJ+1.  In  matrix  form, 

0  =  Rt  +  fil  (3.16) 

where  0  is  the  L  x  1  vector  of  zeros,  1  is  the  Lx  1  vector  of  ones,  and  R  is  the  matrix 

— (//  4-  A)  A  0  0 

H  ~(2/i  4-  A)  A  0 

0  2/i  (3/i  4-  A)  A 

0  ...  0  (L-  2)/i  -((L-l)/i  +  X) 

0  ...  0  0  (L-l)n 

A  similar  analysis  for  finite  transmitters  gives  the  matrix 

-U-V)  0  0 

«  -((M-+J|»)  C*f  -  Z)A  o 

0  ^  (,«%;>)  (Af  -  3)A 

0  ...  0  (L-2)n 

0  ...  0  0 

We  can  solve  Eqn.  3.16  for  •  as 

f  =  -/iR~l 1  (3.19) 

Thus,  having  found  Ps\j,  we  can  find  P5  and  7. 


(M-L  4-1)  A 

(L  -  l)/i  -((a/-2)a)  y 

(3.18) 
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S3  Throughput  Analysis 


Because  it  is  possible  for  a  packet  to  be  successfully  transmitted  even  if  overlapping 
transmissions  occur,  there  are  several  steps  required  in  finding  the  throughput.  For  an 
infinite  population,  throughput  S  is  defined  as 

CO 

S  =  j (Fraction  of  time  spent  transmitting  j  successful  packets)  (3.20) 

j= o 

For  finite  transmitters,  5,  ,  the  throughput  for  user  t  is  the  fraction  of  time  spent  transmit¬ 
ting  successful  packets.  Network  throughput  can  be  found  as  the  sum  of  user  throughputs, 
which  is  simply  A/S,,  since  all  users  are  identical. 

Unfortunately,  the  probability  of  a  packet  being  successful  depends  on  its  length  as 
well  as  the  state  of  the  network  upon  its  arrival.  Thus,  throughput  is  not  simply  7 /fi,  as  it 
would  be  for  fixed  length  packets  (7  is  the  expected  number  of  successful  packets  over  time 
as  in  Section  3. A).  Furthermore,  for  our  definition  of  success,  a  packet  cannot  be  counted 
as  successful  until  it  is  completed,  as  an  error  anywhere  will  cause  the  entire  packet  to  be 
unsuccessful.  Thus,  for  determining  success  or  failure,  it  is  not  sufficient  to  know  only  the 
state  of  the  Markov  chain  at  the  beginning  of  a  transmission.  For  this  reason,  we  follow  the 
approach  used  by  Brazio  and  Tobagi  [BRAZ84J,  and  introduce  an  auxiliary  Markov  chain 
which  traces  the  evolution  of  a  tagged  packet  until  either  an  error  occurs  or  the  packet 
transmission  is  successfully  completed.  We  define  S'  to  be  the  subset  of  the  state  space 
S  consisting  of  all  states  from  which  a  successful  transmission  can  begin.  An  arrival  to  a 
state  in  S'  will  cause  a  transition  to  a  state  in  the  auxiliary  Markov  chain. 

We  define  the  following  parameters: 

Ps\j,t  i3  the  probability  that  a  packet  is  successful  given  that 

i)  it  is  of  length  t, 
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ii)  the  radio  is  not  busy  when  the  packet  arrives, 


iii)  it  finds  the  network  in  state  j  upon  arrival,  and 

iv)  the  protocol  and  channel  state  j  allow  it  to  be  transmitted  (in  the 
channel  load  sense  model  discussed  in  Section  6). 

Tt  is  the  successful  length  of  a  packet,  which  is  0  if  the  packet  is  unsuccessful,  and  is 
the  packet’s  length  if  it  is  successful. 

E(Tg\j)  is  the  expected  successful  length  of  a  packet  given  that  the  network  is  in  state 
j  upon  its  arrival,  and  that  the  packet  was  transmitted. 

Pi, dle  U)  the  probability  that  user  i  is  not  transmitting  given  state  j,  as  in  Section 
3. A. 

E(Tg\j)  is  found  from  Ps\j,r  by  removing  the  condition  on  the  packet  length.  Thus, 

E(r.U)  =  jH  Psy,  r/r(r)  dr  (3.21) 

where  /j(t)  is  the  pdf  of  the  packet  lengths.  Here,  /^(t)  =  ne~^T .  We  note  that  Ps\j 
derived  in  Section  3. A  is 

ps\j  =  JQ  P s|y,r/r(r)dr  (3.22) 

For  an  infinite  population,  throughput  is  given  by 


oc  L—  1 

s  =  £  *«jE{T,\i)  =  £  Airy£(r.|j)  (3.23) 

3= 0  3=0 

As  shown  in  Appendix  A,  for  finite  transmitters,  user  throughput  is  given  by 

A/-1 

5,  =  E  ^3Pi,oLEij)E{T.\j)  (3.24) 

3-0 
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As  in  Section  3. A, 


M  -  i 

piiDLEU)  =  —fif-  for  *  =  !» 2»  •••> M  (3-25) 

So,  network  throughput  5  is 

min(L,A/)— 1 

S  =  AfS{  =  £  \x,iM  -  j)E(T.\j)  (3.26) 

}=o 


3.C  Calculation  of  E(Tt\j) 

The  simple  case  L  =  1,  infinite  transmitter  population,  was  analyzed  by  Ferguson 
[FERG75]  and  by  Bellini  and  Borgonovo  [BELL80].  The  probability  of  a  packet  of  length 
r  being  correct,  PS|0t,  is  the  probability  of  no  new  arrivals  in  time  r,  or 

Ps\o,r  =  e'Ar  (3.27) 


From  equation  3.21, 


E{T,\0)  =  jf°VAr  rue-^dr 

M  =  1/M 

(A+/02  (1  +  ?)2 


which  gives 


5  = 


Ae"* 


j/M 

(1  +  ff)2 


9*-« 

(1  +  <7)2 


The  maximum  throughput  of  0.137  occurs  at  g  —  \/2  -  1. 


(3.28) 


(3.29) 


To  solve  for  E(Tt\j)  for  L  >  1,  we  consider  the  evolution  of  a  packet  that  arrives  and 
is  transmitted.  If  j  >  L,  E(Tt\j)  =  0.  If  j  <  L ,  the  packet  will  begin  to  be  transmitted 
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with  a  non-zero  probability  of  success.  As  the  packet  transmission  progresses,  the  network 
state  changes.  Eventually,  either  an  error  occurs,  in  which  case  the  transmission  fails,  or 
the  packet  is  completed,  and  the  transmission  succeeds.  This  evolution  is  modeled  by  an 
auxiliary  Markov  chain  with  a  state  space  $,,**,  which  is  constructed  from  $  as  follows 
[BRAZ84j. 

First,  we  delete  all  states  for  which  the  probability  of  successful  transmission  is  zero. 
For  an  infinite  population,  this  leaves  states  {0,1,... ,£},  and  for  finite  transmitters,  the 
states  {0, 1,  ...,min(M,  L)}.  If  M  <  L,  all  packets  will  be  successful,  so  E[T,\j)  =  1/p,  the 
average  packet  length.  Thus,  we  will  consider  only  cases  where  L  <  M.  Next,  we  delete  all 
states  corresponding  to  no  transmissions  occurring,  here  the  state  j  =  0.  Finally,  we  add 
two  absorbing  states,  Success  and  Failure.  We  reduce  by  fj,  the  rate  of  all  transitions  from 
states  j  to  states  j  —  1,  2  <  j  <  L,  and  add  transitions  from  every  non-absorbing  state 
to  the  state  Success  at  rate  ft.  A  transition  from  the  state  L  to  the  state  Failure  is  added 
corresponding  to  the  rate  at  which  new  transmissions  are  scheduled  when  in  the  state  L , 
which  is  A  for  an  infinite  population,  and  is  (M  —  L) A  for  finite  transmitters.  Thus,  we 

have  the  auxiliary  Markov  chains  given  in  Figs.  3.3  and  3.4.  If  we  index  the  states  Success 

and  Failure  as  the  last  two  states,  the  transition  rate  matrix  R*  is 

/R  f il 

R*  =  0  0  0  (3.30) 

V  0  0  0  / 

where  R  is  the  LxL  sub-matrix  corresponding  to  transitions  between  the  states  {1,2, ...,  L } 
of  the  auxiliary  Markov  chain,  y?  is  the  lx  L  sub-matrix  corresponding  to  transitions  to  the 
state  Failure,  1  is  the  1  x  L  vector  of  ones,  and  0  is  the  Lx  1  vector  of  zeros.  For  infinite 
and  for  finite  transmitters,  R  is  identical  to  the  corresponding  matrix  given  in  Section 
3. A.  Indeed,  as  we  show  in  Appendix  B,  P$\j  is  simply  the  probability  that  the  auxiliary 
Markov  chain  is  in  state  success  at  time  r  =  oo  given  that  it  was  in  state  j  +  1  at  time 
r  =  0,  which  can  be  found  from  •  =  -/xR_1l  as  shown  earlier.  Furthermore,  the  quantity 
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Figure  3.3.  Auxiliary  Markov  chain  for  the  infinite  population  model. 


IS 


Figure  3.4.  Auxiliary  Markov  chain  for  the  finite  transmitters  model. 


E(Tt\j)  is  the  expected  time  to  reach  state  Success  given  that  the  auxiliary  Markov  chain 
starts  in  state  j  at  time  r  =  0.  This  has  been  shown  [BRAZ84]  to  be 

E{T.\J)  =  [^E-2l]y+,  i  e  S'  (3.31) 

The  derivation  is  given  in  Appendix  B. 


3.D  Packet  Header  Synchronisation  Probability 


Because  this  model  does  not  specifically  account  for  the  receivers,  packet  header  syn¬ 
chronization  affects  the  probability  of  a  packet  being  successful  but  it  does  not  affect  the 

state  space.  We  define  ay  to  be  the  probability  that  the  receiver  captures  the  packet  given 

j  interfering  transmissions.  The  contribution  to  throughput  by  packets  that  found  the 

network  in  state  j  upon  arrival  will  be  reduced  by  the  factor 

otj.  Similarly,  7  and  Ps  will 

be  reduced.  For  an  infinite  population,  we  have 

L- 1 

ps  =  E  *3a3PS  1/ 

3=0 

(3.32) 

L- 1 

I  =  £  Xn3a3PS\j 

3=0 

(3.33) 

S  =  Z  XxjCtjE(Tt\j) 

3=0 

(3.34) 

and  for  finite  transmitters,  we  have 

(3.35) 

A/-1 

7  =  £  Xnj(M  -  j)otjPs\j 

3=0 

(3.36) 

S  =  MZ  X*j{M  -  j)ajE(T.\j) 

(3.37) 

3=0 


Again,  if  L  <  M,  PS\3  and  E{Tt\ J)  are  zero  for  j  >  L. 
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S.E  Results 


In  Figs.  3.5  and  3.6  we  plot  the  throughput  S,  7//i,  and  Ps  for  L  =  10  for  an  infinite 
population,  and  for  finite  transmitters  ( M  =  25).  These  results  are  similar  to  those  found 
by  Pursley  [PURS83].  We  notice  that  when  Ps  drops  below  one,  7/ju  becomes  greater 
than  the  throughput,  but  that  the  curves  are  similar. 

In  Fig.  3.7,  we  plot  normalized  throughput  S/L  vs.  normalized  offered  traffic  g/L 
for  several  values  of  L  for  the  infinite  population  model.  Although  L  is  limited  to  whole 
numbers,  a  continuous  function  is  shown  in  the  following  plots  to  make  them  easier  to 
read.  The  results  are  analogous  to  those  found  by  Musser  and  Daigle  [MUSS82]  for  fixed 
size  packets.  We  notice  that  as  L  increases,  the  maximum  throughput  comes  closer  to  the 
capacity  Z,  but  performance  is  more  sensitive  to  the  offered  traffic  rate. 

4.  Realistic  Channel  Model 

The  performance  of  spread  spectrum  radio  channels  using  error  correction  coding  in  a 
CDMA  environment  has  been  analyzed  by  many  authors.  The  models  derived  either  bound 
or  approximate  the  channel  performance.  Several  difficulties  arise  in  incorporating  these 
models  into  the  Markovian  network  model  described  earlier.  First,  the  detailed  channel 
models  require  a  large  amount  of  computation  for  even  moderate  numbers  of  interfering 
users.  Thus,  numerical  results  are  not  tractable  over  the  entire  range  of  interest.  Second, 
the  performance  of  decoders  is  often  modeled  using  somewhat  loose  bounds.  These  bounds 
are  poor  approximations  to  actual  performance  for  regions  of  high  BER.  However,  these 
high  BER  regions  are  within  the  range  of  interest  for  network  performance.  Finally,  the  use 
of  decoders  introduces  memory  into  the  system,  as  the  output  bit  depends  on  the  received 
signal  over  an  interval  extending  many  bits  into  the  past.  In  this  section,  we  discuss  the 
various  approximations  used  to  overcome  these  difficulties  and  verify  the  validity  of  the 
approximations.  We  then  present  numerical  results  for  several  cases. 
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Throughput  and  gamma/mu 


Figure  3.5.  Throughput,  7 //z,  and  Ps  versus  offered  traffic  g 
for  infinite  population  and  L  =  10. 


Throughput  and  gamma/mu 


Offered  traffic  g 


Figure  3.6.  Throughput,  7//*,  and  F5  versus  offered  traffic  g 
for  finite  transmitters  (M  =  25)  and  L  =  10. 
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Figure  3.7.  Normalized  throughput  S/L  versus  normalized  offered  traffic  g/L 

for  an  infinite  population. 


4.A  Radio  Channel  Model  [TAIF84] 


As  mentioned  previously,  a  large  number  of  detailed  analyses  of  DS-CDMA  channels 
have  been  published.  We  choose  to  incorporate  a  model  which  is  numerically  tractable 
even  for  large  numbers  of  interfering  users.  The  model  was  developed  and  analyzed  by 
Taiple  [TAIP84].  It  relies  on  several  worst  case  assumptions,  and  does  not  require  specific 
characteristics  of  the  set  of  codewords  used  for  the  PN  sequence.  Therefore,  some  of  the 
approximations  used  are  lower  bounds  which  are  believed  to  be  only  moderately  tight  over 
some  of  the  range  of  interest. 

The  channel  model  is  of  a  coherent  BPSK  DS-CDMA  communications  system.  It  does 
not  account  for  fading,  multipath,  or  jamming.  The  received  signal  power  is  identical 
for  all  transmitter-receiver  pairs,  and  the  thermal  noise  level  is  identical  for  all  users,  so 
all  receptions  have  the  same  received  bit  energy  to  noise  energy  E^/Nq.  The  receiver  is 
assumed  to  be  perfectly  synchronized  with  the  transmitted  signal,  which  includes  carrier 
phase,  pseudo-noise  code  chip  timing,  and  bit  timing.  In  addition,  a  worst  case  is  assumed 
for  the  interfering  signals;  all  interfering  signals  are  in  phase  and  chip  aligned  with  the 
desired  signal. 

We  adopt  the  following  terminology.  Bits  refer  to  the  uncoded  information  stream 
generated  by  the  user.  Symbols  refer  to  the  output  of  the  error  correction  encoder.  Thus, 
for  a  rate  1/2  code,  there  will  be  2  symbols  per  bit.  Chips  refer  to  the  pseudo-noise  coded 
information. 

The  general  channel  model  is  as  follows.  Information  to  be  transmitted  is  formed  into 
packets  of  data  bits.  FEC  coding  is  applied  to  these  bits,  resulting  in  binary  symbols.  Each 
symbol  is  multiplied  by  a  pseudo-noise  (PN)  code,  resulting  in  a  number  of  binary  chips. 
We  denote  the  number  of  chips  per  data  bit  as  N.  N  is  the  total  bandwidth  expansion 
due  to  both  FEC  coding  and  PN  spreading.  Finally,  the  chips  are  input  to  a  modulator 
which  phase  modulates  an  RF  carrier.  Thermal  noise  and  a  number  of  interfering  signals 
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are  added  to  the  transmitted  carrier,  resulting  in  a  corrupted  received  signal.  The  receiver 
coherently  demodulates  the  signal  and  despreads  it  by  multiplying  by  the  PN  code  stream. 
It  then  integrates  over  one  symbol,  and  makes  a  hard  decision  on  the  symbol  value.  We 
denote  the  symbol  time  by  T.  The  decoder  then  estimates  the  source’s  data  bit  stream 
from  the  corrupted  symbol  stream  received.  If  any  bit  is  incorrectly  estimated,  we  declare 
the  entire  packet  in  error.  We  again  implicitly  assume  additional  error  detection,  such  as  a 
CRC  code,  but  ignore  the  overhead  required.  The  radio  and  channel  model  is  diagrammed 
in  Fig.  4.1. 

The  mathematical  model  is  a  special  case  of  the  channel  model  given  by  Pursley 
[PURS77].  All  users  are  identical,  so  for  notational  convenience,  we  assume  that  the 
source  is  user  1  and  the  destination  is  user  0.  In  the  general  model,  the  received  signal 
r(f)  is  given  by 


J 

r(0  =  n(t)  +  £  \^Pak{t  ~  rk)bk{t  -  rk)  cos (uct  +  <pk)  (4.1) 

*=i 

where  J  is  the  number  of  signals  simultaneously  arriving  at  user  0  at  time  t}  ak  is  the 
pseudo-noise  code  stream  of  transmitter  k,  bk  is  the  FEC  coded  symbol  stream  of  trans¬ 
mitter  fc,  4>k  is  the  phase  shift  from  transmitter  k  to  user  0,  rk  is  the  time  delay  from 
transmitter  k  to  user  0,  P  is  the  received  signal  power,  and  n(t )  is  the  additive  white 
Gaussian  noise  (WGN). 

In  the  specific  channel  model  examined,  we  assume  that  the  symbols  bk  and  the  code 
chips  ak  are  independent  random  variables  assuming  the  values  +1  and  —1  with  probabil¬ 
ities  1/2  and  1/2.  In  addition,  because  of  the  worst  case  assumption,  rk  is  a  multiple  of 
the  chip  time  T/N  and  <f>k  is  a  multiple  of  2tt. 

Because  of  the  symmetry  of  the  channel,  the  probability  of  error  is  the  same  whether 
the  transmitted  symbol  is  a  +1  or  a  —1.  Thus,  we  will  only  analyze  the  case  when  the 
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Figure  4.1.  Spread  spectrum  radio  model. 
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transmitted  symbol  is  a  +1.  We  define  ck(t)  as  the  product  of  the  received  chip  from  user 
k  with  the  code  stream  of  user  1,  so 

Ck(t)  =  *k(t  ~  n)bk(t  -  Tk)ax(t  -  rk ) 


Because  of  the  assumption  of  chip  alignment,  we  can  unambiguously  define  the  discrete 
chips 


l  =  1  to  N 


The  random  variables  {c*  /}  are  jointly  independent  Bernoulli  trials. 


The  output  of  the  correlation  receiver  is 


■Hi: 


r(0a*(0  COS (u>ct)dt 


(4.2) 


This  reduces  to 

Z  =  V^fl+EE‘*,l|+'W  (4.3) 

V  k—2l=l  ) 

where  Et  =  PT  is  the  energy  per  symbol,  and  M  is  a  Gaussian  random  variable  with 
variance  N0/2.  It  can  be  seen  that  Z  =  P  +  I  +  M ,  where  P  is  the  desired  symbol,  J  is 
the  sum  of  the  interfering  signals,  and  M  is  thermal  noise.  A  symbol  error  occurs  if  Z  <  0. 
Knowing  D  +  I ,  we  can  find  Pr(Z  <  0)  from  standard  communications  theory.  But  P  is 
known  to  be  s/W,  and  I  has  a  known  distribution,  so  we  can  find  the  mean  probability  of 
symbol  error  given  J  simultaneous  transmissions,  Pe,*ymboi{J),  by  averaging  over  all  values 
of  J. 

=  =  0  <i<U-l)N  (4.4, 


P„^UJ)  =  "e'V  ((j  -  l)JV)2-'J-1'".  (i) 
1=0  \  1  / 


(4.5) 
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where 


P(0  =  Q  (sJlEJNo  (1  +  (4.6) 

and  Q( x)  is  the  Marcum  Q  function, 


(4.7) 


4.B  Error  Correction  Coding  Model 

The  performance  of  CDMA  radio  networks  can  often  be  improved  by  introducing 
FEC  coding.  We  analyze  the  network  performance  using  a  rate  1/2  constraint  length  7 
convolutional  code  with  Viterbi  decoding.  The  decoder  performance  analysis  is  discussed 
in  some  detail  in  Appendix  C.  In  this  section,  we  will  briefly  discuss  the  decoder  model  and 
several  limitations.  From  the  analysis  of  the  decoder,  we  can  find  the  so-called  probability 
of  first  error,  Pe,  which  is  the  probability  of  a  bit  being  in  error  given  that  no  earlier  bits 
were  in  error.  Since  we  declare  an  entire  packet  to  be  in  error  if  any  bit  is  in  error,  this  is 
the  appropriate  performance  measure.  Fig.  4.2  shows  Pe  vs.  the  input  symbol  error  rate. 

Taiple  [TAIP84]  shows  that  for  fixed  length  packets  of  length  £.  in  a  slotted  system, 
the  probability  of  a  packet  being  correct  is  lower  bounded  by  (1  —  Pe)1' ■  For  packets  of 
1000  bits,  if  Pe  is  as  high  as  10-3,  there  is  still  at  least  a  37%  chance  that  the  packet  will 
be  correct.  For  Pe  =  5  x  10~3,  this  drops  to  at  least  0.67%.  This  indicates  that  we  should 
strive  for  a  decoder  model  which  is  accurate  up  to  Pe  of  at  least  5  x  10~3.  Unfortunately, 
as  discussed  in  Appendix  C,  the  analysis  uses  a  union  bound  which  becomes  very  loose  for 
Pe  above  about  10~3.  The  rapid  increase  in  Pe  in  the  range  10-3  to  10-2  can  be  seen  in 
Fig.  4.2.  Thus,  the  decoder  model  gives  pessimistic  results  for  Pe  in  this  range. 
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Probability  of  First  Error 


Figure  4.2.  Probability  of  first  error  versus  symbol  error  rate. 
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4.C  Refined  Network  Model 


From  the  channel  model  and  the  FEC  decoder  model,  we  can  find  P^ij),  the  probability 
of  first  error  given  a  symbol  error  rate  of  Pe,iymboi{j)-  Pe{})  vs.  j  is  plotted  in  Fig.  4.3 
for  Efr/No  =  8.0  for  several  values  of  N. 

As  shown  in  Appendix  D,  for  a  packet  which  does  not  visit  the  state  L  +  1  and  spends 
Lj  bits  in  state  jt  the  probability  of  correct  reception  is  approximately  bounded  by 

Pci  fill- PEU)f‘  (4.8) 

J=  1 

We  can  account  for  this  in  the  network  model  by  introducing  a  time  varying  Poisson  process 
that  generates  errors  at  a  rate  ey  while  the  tagged  packet  is  in  state  j.  Specifically,  this 
is  accomplished  by  adding  transitions  from  each  non-absorbing  state  j  in  the  auxiliary 
Markov  chain  to  the  state  Failure  at  at  rate  ey,  as  shown  in  Fig.  4.4.  The  probability 
that  no  errors  occur  while  in  state  j  for  time  ry  is  simply  e~V>.  For  a  small  bit  time  At, 
Tj  m  CjAt  for  some  number  of  bits  £y.  The  probability  that  a  packet  is  successful  is  the 
product  of  the  probabilities  that  no  errors  occur  during  any  state  j . 

L  L 

Pc=  ne~V'  *  ri(e"A'tj)£y  (4-9) 

/= 1  J=1 

If  we  set  1  —  PeU)  =  e~^uJ ,  the  probability  of  a  packet  being  correct  will  be  precisely  the 
bound  given  by  Eqn.  4.8.  Thus, 

=  ^  =  -ln(1--££-°!)  =  -»'■>(>  -  PeU)),  (4.10) 

here  b  =  l/(/iAt)  is  the  average  number  of  bits  per  packet. 

Several  approximations  are  implied  by  this  model. 
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Figure  4.3.  Probability  of  first  error  versus  number  of  simultaneous  transmissions. 
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»)  First,  the  channel  and  decoder  models  yield  an  error  process  which  occurs  at  bit 
boundaries.  There  are  two  approaches  to  incorporating  this  into  the  network  model.  One 
approach  is  to  change  the  network  model  to  a  discrete  time  model,  with  arrivals  and 
departures  occurring  at  bit  boundaries.  For  bit  times  which  are  short  in  comparison  to 
the  average  holding  times  in  the  states,  this  model  will  be  very  similar  to  the  continuous 
time  model,  except  that  there  will  now  be  transitions  occurring  between  non-neighboring 
states.  These  transitions  will  have  very  low  probabilities  for  short  bit  times,  as  the  arrival 
or  departure  rate  per  bit  is  low.  We  chose  the  alternative  approach,  which  is  to  approximate 
the  error  process  as  a  continuous  time  process,  in  which  errors  are  generated  by  a  time 
varying  Poisson  process.  This  yields  a  simpler  network  model  and  will  give  accurate  results 
for  average  packet  lengths  of  1000  bits. 

«i)  Second,  the  network  model  requires  a  cut-off  L  such  that  errors  occur  with  probabil¬ 
ity  one  for  j  >  L.  No  such  cut-off  exists  in  the  physical  channel,  since  the  probability  of  bit 
error  is  never  greater  than  1/2,  so  there  is  always  a  non-zero  probability  of  correct  recep¬ 
tion.  Introducing  the  limit  L  underestimates  the  throughput,  as  some  successful  packets 
will  be  counted  as  unsuccessful.  Nevertheless,  if  L  is  chosen  such  that  (1  —  Pe{L  +  1))^ 
is  small  for  H  equal  to  the  mean  holding  time  in  state  L  4-  1,  the  probability  of  a  packet 
successfully  surviving  a  visit  to  state  L  +  1  will  be  small. 

The  choice  of  L  is  simple  in  the  case  of  the  coded  channel  because  of  the  decoder 
model  used.  As  explained  in  Appendix  C,  the  bound  for  the  decoder  model  is  not  known 
to  be  accurate  for  Pe  above  10~2.  For  this  reason,  the  model  effectively  truncates  the 
performance  in  this  range,  and  sets  Pe  =  1  when  the  bound  indicates  Pe  >  10~2.  Thus, 
L  is  chosen  to  be  the  largest  j  such  that  PeU)  <  10-2. 

m)  Another  approximation  implied  by  the  model  is  that  the  addition  of  an  error  process 
does  not  affect  the  existing  transition  rates  for  the  auxiliary  Markov  chain.  In  reality,  the 
existing  transition  rates  leaving  state  j  will  be  reduced  by  1  -  Pe{])  per  bit.  However,  for 
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Pe{))  as  high  as  10  2,  this  will  only  be  a  1%  reduction.  Thus,  this  approximation  has  a 
negligible  effect  on  throughput  for  the  coded  channel  model. 

iv)  The  final  assumption  made  in  incorporating  a  realistic  channel  model  into  the 
network  model  is  that  the  channel  behavior  is  memoryless.  It  can  be  shown  that  for 
some  of  the  cases  considered,  the  decoder  performance  depends  on  past  states  most  of  the 
time.  Nevertheless,  the  approximate  bounds  described  below  show  that  the  variation  in 
throughput  due  to  this  assumption  is  small.  Thus,  even  though  the  model  does  not  exactly 
parallel  the  physical  system,  the  results  closely  match  the  ideal  performance  of  the  system. 

The  decoded  memory  typically  extends  30  bits  into  the  past  [CLAR81].  Thus,  Pe  at 
time  r  depends  not  only  upon  the  current  state,  but  also  upon  the  states  visited  during 
the  last  30  bits,  or  60  symbols.  The  arrival  and  departure  processes  are  both  Poisson,  so 
the  holding  time  in  a  state  j  is  exponentially  distributed,  with  rate  A  +  j'/x  for  the  infinite 
population  model,  or  (M  -  j)A  +  j/x  for  the  finite  transmitters  model.  The  mean  holding 
time  in  state  j  is  thus 


__I_  =  lliL 
*  +;>  g  +  j 

1 _ _  V/* _ 

(M-j)X+jn  (M  -  j)g  +  j 

For  \/n  =  1000  bits,  the  average  holding  time  is 


infinite  population 
finite  transmitters 


(4.11) 


1000  1000 

- -  QJ*  - - - - - 

g  +  j  (M  -  j)g  +  j 

For  the  coded  channel,  with  Eb/N0  =  8.0,  N  =  64  chips  per  bit,  the  cut-off  L  =  7,  and  the 
maximizing  values  of  offered  traffic  g  are  3.01  (infinite  population),  0.04  (80  transmitters), 
and  0.53  (10  transmitters).  It  can  be  seen  that  even  for  the  state  j  =  L,  the  average 
holding  times  are  1000/(7  +  3.01),  1000/(7 +  80  0.04),  and  1000/(7+ 10  0.53),  all  of  which 
are  about  three  times  the  decoder  memory  of  30  bits.  In  these  cases,  we  can  ignore  the 
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transient  behavior  which  occurs  during  the  30  bits  following  a  state  change.  Unfortunately, 
for  larger  N ,  the  mean  holding  time  in  some  states  is  as  low  as  10  or  20  bits.  For  these 
cases,  the  decoder  output  is  very  likely  to  depend  upon  the  previous  state  for  a  large 
fraction  of  the  time  spent  in  these  states. 

Nevertheless,  even  for  these  few  extreme  cases,  the  probability  that  the  number  of 
interfering  transmissions  j  varied  widely  during  the  last  30  bits  is  not  large.  For  offered 
traffic  g  in  the  range  which  gives  maximum  throughput,  most  of  the  steady  state  probability 
is  concentrated  in  the  states  with  intermediate  values  of  Pe{])-  We  find  that  PeU)  varies 
slowly  with  j  over  these  intermediate  ranges  (see  Fig.  4.3).  This  is  especially  true  for  very 
large  numbers  of  chips  per  bit,  which  are  the  cases  where  the  memory  may  be  significant. 
Thus,  for  those  cases  were  the  decoder  memory  cannot  be  ignored  as  a  second  order 
transient  effect,  the  difference  between  PeU),  PeU  +  1)  and  Pe{ j  —  1),  and  even  Pe[ j  +  2) 
and  Pe(j  -  2)  is  small. 

To  verify  this  approximation,  we  calculated  approximate  bounds  by  using  first  Pe{j+3) 
and  then  Pe{ j  —  3)  as  the  error  rates  from  state  j.  This  was  done  for  the  infinite  population 
model,  with  Eb/No  =  8.0  and  with  256  chips  per  bit.  For  this  rasp  t^p  cut-off  L  is  26. 
For  offered  traffic  g  of  26.0,  the  transitions  from  j  to  j  +  1  occur  at  a  rate  gn  =  26.0/1000 
per  bit,  and  the  transitions  from  j  to  j  —  1  occur  at  a  rate  jfi,  which  is  upper  bounded  by 
Ln  =  26/1000  per  bit.  The  probability  of  going  from  state  j  to  a  state  j'  >  j  +  6  in  30 
bits  is  less  than  the  probability  of  more  than  6  arrivals  from  the  Poisson  arrival  process 
(which  is  rate  0.026  per  bit)  in  30  bits,  which  is  simply 


l  e-°-78(0.78 y 

t=0 


t! 


(4.12) 


For  6  =  3,  this  comes  to  0.0083.  This  probability  is  even  smaller  for  offered  traffic  g  less 
than  26.0.  Similarly,  the  probability  of  going  from  state  j  to  a  state  j '  <  j  —  6  in  30  bits 
is  less  than  0.0083.  Thus,  when  the  current  network  state  is  j,  with  probability  greater 
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than  99%,  the  minimum  symbol  error  rate  encountered  by  any  of  the  60  symbols  in  the 
decoder  memory  will  be  no  less  than  PeU  —  3)  and  the  maximum  symbol  error  rate  will 
be  no  greater  than  Pe{ j  +  3).  Thus,  substituting  Pe(j  ±  3)  for  the  symbol  error  rate  of 
state  j  will  give  approximate  bounds  to  the  performance  of  a  decoder  with  a  memory  of 
30  bits. 

These  bounds  are  fairly  loose,  since  for  many  states,  the  variation  will  be  less  than 
±2  or  ±1  states  with  a  very  high  probability.  Also,  even  when  the  states  j  +  3  or  j  —  3 
are  visited,  many  of  the  last  60  symbols  will  have  symbol  error  probabilities  closer  to 
Pe.»ymboi(j)  than  Pe,gymboi(j  ±  3).  Even  so,  as  can  be  seen  in  Fig.  4.5,  the  maximum  values 
of  the  throughput  for  the  models  which  use  Pe{]  ±  3)  are  within  20%  of  the  maximum 
throughput  for  the  model  which  uses  PeU)-  Therefore,  the  approximate  model  which 
ignores  the  memory  of  the  decoder  yields  meaningful  results  even  for  extreme  cases  where 
holding  times  are  on  the  order  of  30  bits. 

4.D  Standard  ALOHA  Channel 


For  comparison,  we  derive  the  throughput  of  the  uncoded,  unspread  radio  channel 
with  thermal  noise.  In  this  case,  N  =  1  chip  per  bit,  and  also  L  —  1.  The  analysis  is 
similar  to  the  case  with  no  thermal  noise,  which  was  presented  in  Section  3.C.  However,  in 
addition  to  new  arrivals,  there  are  also  failures  due  to  an  error  arriving  before  the  packet 
is  completed.  Thus, 

^5|0,r  =  e_(A+(,)r  (4.13) 


£(T.|0) 


_ 1/M _ 

(1  +  <7  +  e  i/m)2 


(4.14) 


S  -- 


ge  9 

(1  +  9)’ 


(4.15) 
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Throughput 


Figure  4.5.  Throughput  for  bounds  on  the  probability  of  symbol  error. 
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For  an  uncoded  channel  with  Eb/N0  =  8.0,  1000  bits  per  average  packet,  ei//i  =  1.15x  10~5, 
which  is  much  less  than  1,  so  the  extra  term  is  negligible.  Thus,  the  thermal  noise  has  little 
effect  on  throughput,  so  the  maximum  throughput  is  again  0.137,  achieved  at  g  =  >/2  —  1. 

4.E  Results 

For  a  fair  comparison  of  throughput,  we  state  most  results  in  terms  of  S/N,  which 
is  throughput  divided  by  N,  the  number  of  chips  per  information  bit.  S/N  indicates  the 
throughput  per  unit  bandwidth,  where  a  bandwidth  of  one  is  required  by  an  uncoded 
unspread  signal. 

Fig.  4.6  shows  the  normalized  throughput  S/N  vs.  the  normalized  offered  traffic  g/N, 
for  the  coded  channel  with  E^/No  =  8.0,  for  the  infinite  population  model.  Results  are 
given  for  several  values  of  N.  This  has  a  shape  similar  to  Fig.  3.7. 

The  following  curves  show  the  effect  of  various  parameters  on  the  maximum  throughput 
5mlx,  achieved  by  maximizing  over  offered  traffic  g.  For  convenience,  we  refer  to  this 
maximum  as  throughput  or  normalized  throughput,  without  specifically  indicating  the 
maximization  over  g. 

Fig.  4.7  shows  the  throughput  vs.  N  for  several  values  of  E^/No,  for  the  infinite 
population  model.  The  corresponding  values  of  S/N  are  plotted  in  Fig.  4.8.  It  can  be 
seen  that  the  performance  becomes  limited  by  multi-user  interference  for  E^/Nq  greater 
than  about  20.  Also,  the  normalized  throughput  increases  monotonically  with  N .  The 
maximum  normalized  throughput  for  the  uncoded  ALOHA  channel  ( N  =  1)  is  0.137. 
Thus,  over  all  of  the  ranges  of  E(,/No  and  N  examined,  the  normalized  capacity  of  the 
spread  spectrum  channel  with  FEC  coding  is  less  than  that  of  a  non-spread  channel  with 
no  coding. 

We  plot  S/N  vs.  Ef,/No  for  the  coded  and  the  uncoded  channels  in  Fig.  4.9,  for  256 
chips  per  bit.  The  improvement  due  to  FEC  coding  is  significant.  Again,  we  notice  that 
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Figure  4.6.  Normalized  throughput  S/N  versus  normalized  offered  traffic  g/N 

for  an  infinite  population. 
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Normalized  throughput  S/N  versus  chips  per  bit  N 
for  an  infinite  population. 
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Normalized  Throughput 


Figure  4.9.  Normalized  throughput  S/N  versus  Eb/No 
for  an  infinite  population. 
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the  curves  flatten  out  for  large  Eb/No ,  as  the  multi-user  interference  begins  to  dominate 
the  thermal  noise. 

Fig.  4.10  shows  5  vs.  N  for  the  infinite  population  case  and  for  several  values  of  Af 
in  the  finite  transmitters  case,  for  Eb/No  =  8.0.  As  N  increases,  the  channel  reaches  a 
point  where  Pe{ Af)  becomes  small.  At  this  point,  even  with  all  M  users  transmitting, 
there  are  few  unsuccessful  packets.  Thus,  increasing  N  further  results  in  little  increase  in 
maximum  throughput.  This  can  be  seen  for  the  M  =  10  and  the  M  =  20  curves.  Because 
of  the  additional  bandwidth  required,  increasing  N  can  actually  cause  a  decrease  in  the 
normalized  throughput.  This  is  seen  in  Fig.  4.11.  We  note  that  the  peak  value  of  S/N 
for  M  =  10  and  M  =  20  is  less  than  half  the  throughput  for  the  uncoded  N  =  1  channel. 

In  Fig.  4.12,  we  plot  the  network  throughput  and  also  the  individual  user  throughput 
vs.  the  number  of  users  for  the  finite  transmitters  model,  Eb/No  =  8.0  and  N  =  128  chips 
per  bit.  It  is  interesting  to  find  that  there  is  an  optimum  value  of  M  which  maximizes 
network  throughput.  We  will  see  in  Section  5  that  this  is  no  longer  true  when  the  effect 
of  half-duplex  radios  is  taken  into  account. 

5.  Half-Duplex  Radios 

6.  A  Analysis 

We  next  consider  a  finite  user  network  where  both  the  sources  and  the  destinations  of 
packets  are  among  the  M  users.  Several  modifications  are  needed.  First,  we  must  account 
for  the  case  when  the  source  is  busy  receiving  a  packet,  and  therefore  will  not  transmit  a 
packet  which  is  scheduled.  Second,  we  must  account  for  the  case  when  the  destination  is 
busy  receiving  or  transmitting,  and  will  not  receive  the  packet.  Finally,  we  must  treat  the 
effect  on  state  transitions  of  an  idle  receiver  not  capturing  a  packet. 

The  state  variable  J C(r)  is  now  defined  as  both  the  number  of  active  transmitters  and 
the  number  of  receivers  which  have  captured  a  packet.  This  gives  a  state  space  S  =  {(<,  r)  : 
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Throughput 


Figure  4.10.  Throughput  S  versus  chips  per  bit  N 
for  an  infinite  population  and  for  finite  transmitters. 
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Figure  4.12.  Throughput  versus  user  population  M 
for  finite  transmitters. 
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0<t<M,0<r<t ,  and  t  +  r  <  M }.  For  the  basic  model,  we  observed  that  successful 
transmissions  could  only  begin  when  the  state  was  in  a  subset  $'  of  the  state  space  $} 
corresponding  to  j  <  L.  In  the  model  with  finite  receivers,  we  find  that  this  subset  isjiow 
S'  =  {(*>  r)  :  0  <  t  <  M  —  l,  t  <  L,  0  <  r  <  t,  and  t  +  r  <  M  —  1}.  The  non-absorbing 
states  of  the  auxiliary  Markov  chain  consist  of  the  subset  {(f,  r)  :  (t  —  l,r  —  1)  G  S'}. 

We  define  the  following  variables: 

Ps\(t,r)  is  the  probability  that  a  transmitted  packet  is  successful  given  that  the  network 
was  in  state  (f ,  r)  upon  its  arrival  and  given  that  it  was  captured. 

Psi(t,r)  is  the  probability  that  t^e  source  radio  is  idle  given  state  (t,r) 

Pj}l{tt  r)  is  the  probability  that  the  destination  radio  is  idle  given  state  (<,  r)  and  given 
that  the  source  is  idle 

at  is  the  probability  that  the  destination  radio  synchronizes  to  the  packet  header  given 
that  it  is  idle  and  given  that  there  are  t  interfering  transmissions. 

rj  is  the  steady  state  distribution  of  Jf(r) 

n rj  is  the  steady  state  probability  of  the  network  being  in  state  ( t ,  r)  at  time  tJ"  given 
that  a  packet  transmission  begins  at  time  ro 

E(T,\(t}  r))  is  the  expected  successful  length  as  in  the  basic  model. 

Pe{ 0  is  the  probability  of  first  error  for  t  simultaneous  transmissions. 

Because  of  identical  users,  for  any  user  i,  we  find 

„  .  .  M  —  t  —  r  , 

Psi(ttr)  =  - — M  (5.1) 

PDl(t,  r)  =  1  -  (5.2) 
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since  the  destination  is  uniformly  distributed  over  the  M  —  1  other  radios,  t  +  r  of  which 
are  busy. 

In  state  {t,  r),  transmissions  are  scheduled  at  rate  (M  —  t  —  r)A.  A  packet  has  a 
probability  atP£>i(t,  r)  of  being  captured,  which  will  cause  a  transition  to  state  (<+1,  r-fl). 
Otherwise,  the  packet  will  not  be  captured,  so  it  will  cause  a  transition  to  state  (f  +  1,  r). 
Of  the  t  packets  being  transmitted,  r  have  captured  a  receiver  and  t  —  r  have  not  captured 
a  receiver.  Thus,  due  to  packets  transmissions  completing,  there  is  a  transition  to  state 
(t  -  1,  r  —  1)  at  rate  r/x,  and  to  state  (t  —  1,  r)  at  rate  ( t  —  r)/x.  The  state  transitions  for  the 
general  state  (f,r)  are  diagrammed  in  Fig.  5.1.  The  forms  of  the  Markov  chain  X{r)  and 
of  the  auxiliary  Markov  chain  are  given  in  Figs.  5.2  and  5.3.  The  steady  state  probabilities 
7T(,  r)  of  the  Markov  chain  X(t)  are  not  easily  found  analytically.  However,  they  can  be 
found  from  the  conservation  equations 

Qir  =  0  (5.3) 

where  Q  is  the  transition  rate  matrix,  and  from 

£  =  1 
(t,r)eS 

Similar  to  the  result  in  Section  3. A,  we  find 

PS  =  £  *{<1r)ar</>^(e>r)PS|(f,r) 

(l.r)eS' 

E  Psi{t,  r)wit  r)atPDI(t,  r)Ps i(t  r) 

_  (<,r)6S'  _ _ _ 

E  Psi{t,  »•)«■(<, r) 

MeS' 

As  we  show  in  Appendix  E, 

Si  =  £  A7r(tir)P5/(<,r)PjD/(t,r)atE(T,|(t,r))  (5.6) 

(t.r)eS' 


(5.4) 
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Figure  5.2.  State-transition  diagram  for  the  half-duplex  model. 
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Since  S  =  MS,,  we  find 


•*  =  £  A^,,r)(M-l-r)(l-±fT)0,£(T.|((,r))  (5.7) 

Psj(t,r)  and  f;(r,|(t,r))  are  again  found  from  the  auxiliary  Markov  chain  transition  rates. 
The  Poisson  error  generation  process  described  in  Section  4  is  included  in  this  auxiliary 
Markov  chain. 

(.B  Results 

Throughput  vs.  N  is  plotted  for  several  values  of  M  in  Fig.  5.4.  We  find  that  the 
throughput  is  much  lower  than  for  the  same  values  of  M  in  the  basic  model.  This  is 
expected,  as  the  maximum  possible  throughput  is  [M/2\.  The  jump  at  the  point  where 
the  curve  flattens  out  occurs  when  L  =  M  —  1.  This  jump  is  due  to  the  truncation  of 
the  decoder  performance  model  for  Pe  >  10~2.  Below  the  jump,  the  decoder  model  uses 
the  loose  bound  Pe  —  1.0  for  some  states,  while  above  the  jump,  the  model  uses  a  tighter 
bound  for  all  states. 

The  maximum  throughput  flattens  out  around  the  point  where  the  channel  can  support 
M/2  simultaneous  transmissions  with  high  probability  of  success.  As  t,  the  number  of  users 
transmitting,  increases  above  M/2,  the  number  of  users  that  can  be  receiving  is  limited 
to  M  —  t,  so  the  probability  that  the  destination  is  idle  becomes  very  small.  Thus,  we 
expect  that  the  optimum  offered  traffic  is  such  that  t  <  [M/2J  for  most  of  the  time,  so 
increasing  N  to  support  more  than  [M/2J  transmissions  should  not  have  much  impact  on 
maximum  throughput.  In  Fig.  5.5,  we  plot  network  throughput  and  also  user  throughput 
vs.  M  for  Eb/No  =  8.0,  N  =  128.  We  find  that,  contrary  to  the  earlier  results  of  Section 
4,  the  network  throughput  for  the  half-duplex  model  increases  monotonically  with  M. 
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Figure  5.4.  Throughput  S  versus  chips  per  bit  N 
for  the  half-duplex  model. 
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Figure  5.5.  Throughput  S  versus  user  population  M 
for  the  half-duplex  model. 


6.  Channel  Load  Sensing 


Many  authors  have  shown  that  for  ALOHA-like  systems,  the  network  performance 
can  be  improved  by  sensing  the  channel  and  blocking  transmissions  when  the  channel 
is  busy.  In  a  CDMA  network  where  multiple  users  can  transmit  simultaneously,  similar 
improvement  is  possible  if  the  radios  can  sense  the  number  of  radios  transmitting  and 
block  transmission  when  the  channel  is  heavily  loaded.  In  addition  to  combatting  multiple 
user  interference,  this  technique  of  channel  load  sensing  may  also  increase  the  robustness 
in  the  presence  of  jamming,  by  causing  the  channel  capacity  to  degrade  gracefully  as  the 
jamming  signal  power  increases. 

6.A  Infinite  Population 

We  use  K  to  denote  the  number  of  sensed  transmissions  at  which  new  transmissions 
are  blocked.  For  K  >  L,  the  state  space  S  is  just  a  truncation  of  the  basic  model  state 
space.  This  has  the  same  state  transition  rates  as  the  state  space  of  the  M/M/K/K  queue, 
which  is  a  /f  server  loss  system  [KLEI75].  The  auxiliary  Markov  chain  is  the  same  as  in 
the  basic  model;  so  for  this  range  of  K  >  L,  the  only  difference  from  the  basic  model  is  in 
the  initial  state  probabilities  {fly},  now  given  by 

_  (a/m)j7j  ? 

nJ  k 

£(A /#«)*/*» 

k= o 

For  very  small  K,  almost  all  transmissions  will  be  successful,  because  Pe{j )  will  be 
small  for  all  ;  6  5.  In  the  limit  as  g  — ►  oo,  as  soon  as  the  end  of  a  transmission  is  sensed 
a  new  packet  transmission  will  begin.  In  the  case  of  zero  propagation  delay,  this  will 
result  in  there  always  being  exactly  K  users  transmitting.  The  probability  of  success  is 
the  probability  that  the  packet  is  completed  before  an  error  occurs.  Thus,  we  find 

PstK-l,,  =  «“*'  (6-2) 

(3 


(6.1) 


1 


(6.3) 


E{Ts\K-i)  = 


(i  +  c*7a02 


so 


Sm&3 r  — 


K 


(1  +  £  KtvV 


(6.4) 


6.B  Half-Duplex  Model 

i)  K  >  M 

Because  there  will  never  be  any  more  than  Af  transmissions  occurring,  a  channel  load 
sense  threshold  K  >  M  will  give  the  same  performance  as  no  channel  load  sensing. 

ii)  K  <  M  and  e^/n  non-negligible 

In  this  case,  the  state  space  is  truncated  at  t  =  K,  while  the  auxiliary  Markov  chain 
remains  the  same.  Consequently,  the  only  difference  from  the  case  studied  in  Section  5  is 
that  the  initial  probability,  distribution  {ir^ r)}  changes.  Thus,  the  solution  can  be  found 
using  the  equations  derived  in  Section  5. 

iii)  K  <  M  and  e^/fi  «  1 

For  this  range  of  K,  almost  all  captured  packets  will  be  successful.  However,  packets 
may  not  be  captured,  due  to  the  destination  radio  being  busy.  Even  so,  we  note  that  when 
K  is  small  compared  to  M,  the  probability  of  no  capture  is  very  small,  so  the  throughput 
is  maximized  by  a  very  large  g.  Unfortunately,  the  numerical  solution  used  for  finding  the 
throughput  becomes  unstable  as  g  becomes  large.  Therefore,  we  use  an  approximation 
of  the  Markov  chain  which  is  valid  for  such  values  of  g.  For  very  large  g,  we  expect  the 
number  of  transmissions  to  be  K,  since  as  soon  as  a  transmission  ends  a  new  scheduling 
point  will  be  generated  and  transmission  of  a  new  packet  will  begin.  However,  in  order  to 
find  the  throughput,  we  must  find  the  distribution  of  the  number  of  captured  receivers. 
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In  this  section,  we  first  justify  that  the  probability  of  there  being  fewer  than  K  active 
transmitters  is  0(l/g).  We  then  use  an  approximation  of  the  Markov  chain  which  contains 
only  the  states  {(<,  r)  :  t  =  K  —  1  or  t  =  K)  to  find  the  distribution  on  the  number  of 
captured  receivers.  E(Ts\[K  —  l,r)  is  found  from  Eqn.  3.21,  to  give 


K 


s  =  E 


nr, 


(*» 


r=o  (1  +fAr/^)2 

We  denote  the  marginal  probability  of  there  being  t  users  transmitting  by  nt,  so 


*i  =  E  *(i,r)  £  =  min(l,  M-t) 

r=0 


(6.5) 


(6.6) 


We  denote  the  expected  number  of  users  receiving  given  t  users  transmitting  by  E(R\t), 


so 


£(fi|()  =  E’-’r^  f  =  ”>in(f,  M-t) 

r=0 


(6.7) 


We  can  find  jrt  from  the  balance  equations  resulting  from  a  cut  between  states  ( t ,  r)  and 
states  ( t  +  1,  r)  (see  Fig.  5.2  ). 


so 


t  t+l 

£(M  -  *  -  =  £(*  +  l)M*(t  +  l,r) 

r=0  r=0 

( M  -  f)Aft  -  =  (t  + 


A  _  (M  -  (t  -  1)  -  E[R\t  -  1))A 

=  — — - — - irt_i 


(6.8) 


(6.9) 


If  we  define  f?(/|f)  as  the  expected  number  of  idle  radios  given  t  users  transmitting,  we 
have 


E{I\t)  =  E(M  -  t  -  R\t) 
=  M  -  t  -  E{R\t) 


(6.10) 
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so 


£(/|t-l)A. 

— Jj,  *•-' 

W  Wn)E(I\u) 

n  -^T.r  1,0 


(6.11) 


For  t  <  M/2,  there  will  always  be  at  least  one  idle  radio,  so  1  <  E(I\t)  <  M  —  t. 
Therefore, 

(\/n)E(I\t) 


t  +  1 


is  0(g)  for  t  <  M/2, 


(6.12) 


which  implies  that  for  K  <  M/2,  *k-\  is  0  and  *K- 2  is  0  (<7~2)*  Thus,  to  first 

order,  for  large  g,  we  can  consider  the  approximation  to  the  Markov  chain  which  is  shown 
in  Figs.  6.1  and  6.2. 


The  two  cuts  indicated  in  Fig.  6.2  give  the  following  local  balance  equations. 

(■ K  -  r)p*(K,r\  =  Ar  (l  -  aK-i  (l  -  ^ly1)) 

A r  =  (K-  r)fiir{Ktr)  +  (r  + 

Where  Ar  is  the  rate  of  probability  flowing  out  of  state  ( K  —  l,r),  so 


(6.13) 

(6.14) 


Ar  =  [M  -  K  -  r  +  l)Airjjt_1  r) 


(6.15) 


From  these  equations,  we  find 


(r  +  l)t*i r(K,,+1)  = 


(K  ~  r)(iiT(K,r)  (K  w 

- {K  -  r)mK-] 


(6.16) 


so 


*{K,r+ 1)  = 


K-r 


r+i  V1  (i  -  *M=r) 


1  n 


(K,t) 


0  <  r  <  K-  1 


(6.17) 


S6 


Figure  6.1.  State-transition  diagram  for  the  half-duplex 
channel  load  sense  model. 
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Figure  6.2.  Local  state-transition-rate  diagram  for  the  half-duplex 
channel  load  sense  model. 


From  these  recursion  equations,  we  can  find  the  probabilities  K(K,r)  iQ  terms  of  K(k,o)>  and 
then  normalize  by  =  1- 

For  the  case  of  K  small  compared  to  M,  and  e^//i  <<  1,  we  have  found  an  approx¬ 
imation  for  the  Markov  chain  which  gives  the  distribution  7r(A»,  from  which  we  find  the 
throughput  S.  For  larger  K,  the  maximum  throughput  is  achieved  at  values  of  g  for  which 
the  approximation  is  no  longer  valid.  However,  for  g  in  this  range,  we  can  use  the  numerical 
solution  developed  in  Section  5. 

6.C  Propagation  Delay 

It  is  well  known  that  carrier  sensing  systems  degrade  as  the  propagation  delay  becomes 
large  relative  to  the  packet  lengths.  To  investigate  this  effect  for  the  CDMA  model,  we 
introduce  an  exponentially  distributed  holding  time  at  the  beginning  of  each  packet.  The 
decision  on  whether  or  not  to  transmit  is  again  based  on  the  number  of  radios  which  can 
be  sensed  transmitting,  but  this  decision  does  not  account  for  users  in  the  holding  state. 
This  model  approximates  the  behavior  due  to  propagation  delay. 

From  the  standpoint  of  a  receiver,  all  other  users  have  a  state  model  which  is  dia¬ 
grammed  in  Fig.  6.3(a).  During  the  holding  state,  the  receiver  cannot  detect  the  transmis¬ 
sion.  From  the  standpoint  of  a  transmitter,  the  holding  time  and  the  transmission  time 
are  reversed.  Thus,  a  packet  transmission  can  be  illustrated  as  in  Fig.  6.3(b).  The  model 
has  the  peculiarity  that  a  radio  is  busy  both  during  the  transmission  time  and  during  the 
holding  time,  so  it  is  as  if  a  radio  contiuu-  to  be  busy  from  the  time  a  packet  transmission 
is  completed  until  the  time  the  packet  finishes  propagating  to  the  receiver.  This  is  shown 
as  the  shaded  region  in  Fig.  6.3(b).  We  will  see  that  this  extra  busy  period  causes  the 
throughput  to  be  reduced  even  for  no  channel  load  sensing. 
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Figure  6.3.  Exponential  holding  time  model  for  non-zero  propagation  delay. 
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6.D  Delay  Model 


We  introduce  the  parameter  i/,  where  \jv  is  the  average  holding  time.  We  also  define 
h  as  v/n,  so  \/h  is  the  average  holding  time  normalized  by  the  mean  packet  length. 
This  parameter  l/h  is  roughly  equivalent  to  the  standard  delay  measure  a,  which  is  the 
propagation  delay  noimalized  by  the  packet  length.  A  radio  can  be  either  idle,  holding, 
transmitting,  or  (for  the  half-duplex  model)  receiving,  as  we  indicate  in  Fig  6.3(a). 

We  incorporate  the  holding  time  into  the  Markov  chains  as  follows.  Let  w  be  the 
number  of  transmitters  in  a  holding  state.  In  the  model  which  does  not  account  for 
receivers,  the  states  are  indexed  by  w  and  t.  The  state  space  is  5  =  {(to,f)  :  0  <  w  < 
M,  0  <  t  <  M  —  w},  and  the  subset  from  which  successful  transmissions  can  begin  is  S'  = 
{( w,t )  :  1  <  w  <  M,  0  <  t  <  M—w,  and  t  <  L}.  For  the  state  (u/,  t),  t  <  A”,  transmissions 
are  scheduled  at  rate  ( M  —  w  -  <)A,  causing  transitions  to  state  ( w  +  l,f).  Users  change 
from  the  holding  state  to  the  transmitting  state  at  rate  w,  causing  transitions  to  state 
(to  — l,  f  4-1).  Transmissions  are  completed  at  rate  tfi,  causing  transitions  to  state  (w,t  —  1). 
The  state  transitions  for  the  general  state  (u>,  t )  are  diagrammed  in  Fig.  6.4.  For  t  >  K, 
there  are  no  transitions  due  to  new  packet  transmissions,  but  the  other  transitions  are  not 
affected. 

For  the  half-duplex  model,  we  have  a  3-dimensional  state  space  indexed  by  w,  t,  and 

r\ 

0  <  xv  <  M 

S  =  •  (to,  t,r)  :  0  <  t  <  M  -  w 

0  <  r  <  M  —  w  ~  t,  and  r  <  t  . 

1  <  w  <  M 

S'  =  •  (to,  t,  r)  :  0  <  t  <  M  —  w,  and  t  <  L 

.  0<r<M-w-t,  and  r  <  t  . 
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Figure  6.4.  Local  state-transition-rate  diagram 
for  finite  transmitters,  non-zero  propagation  delay. 


Once  again,  channel  load  sensing  causes  the  deletion  of  transitions  due  to  new  packet 
transmissions  for  t  >  K. 


The  derivations  of  user  throughput  5,  are  given  in  Appendix  F.  For  the  model  which 
does  not  account  for  receivers,  we  find 

S,  =  E  v*(tc<t)pH{u>,t)°'tE{T'\lw,t))  (6.18) 

(tr.OeS' 

where  Pjj( to,  t)  is  the  probability  that  user  i  is  in  the  holding  state  given  state  (to,  f),  which 
is  simply  w/M.  So  network  throughput  is 


S=  Y,  ur{Wtt)watE[T9  |(u?,t))  (6.19) 

Again,  E[TK\{w,  <))  is  found  from  n E  21  for  t  <  L.  For  the  half-duplex  model,  we  find 

St=  L'*(w,t,r)Pn{™,t,r)PDI{rv,t,r)alE(T'\{w,t,r))  (6.20) 

(ir,<.r)GS' 

where  P[j{w,t,r)  is  the  probability  that  user  i  is  in  the  holding  state  given  state  ( w,t,r ), 
which  is  simply  rv/M,  and  Fp/(iy,t,r)  is  the  probability  that  the  destination  is  idle  given 
state  (tu,f,r)  and  given  that  i  is  holding,  or 


PDl{w,t,r)  =  1  - 


to  -M  +  r  -  1 

M-  1 


(6.21) 


Thus, 


S  =  E 

(tr.f.r)eS' 


0- 


w  +  t  +  r  —  1 
M  -  1 


)a(F(T,|(«;,t,r))  (6.22) 
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6.E  Results 


For  an  infinite  population,  N  =  256,  and  Ei/Nq  =  8.0,  we  plot  S/N  vs.  g  for  various 
values  of  K,  and  also  for  no  channel  load  sensing  in  Fig.  6.5.  For  these  parameters,  the 
cut-off  L  is  26.  It  is  seen  that  for  small  K ,  the  normalized  throughput  increases  towards 
K/N  as  g  increases.  Also,  the  maximum  throughput  for  K  =  25  is  almost  equal  to  the 
maximum  for  no  channel  load  sensing.  We  then  xdot  the  throughput  5  vs.  the  channel 
load  sense  point  K  for  an  infinite  population  in  Fig.  6.6,  varying  N  from  64  to  1024.  For 
N= 64,  channel  load  sensing  can  increase  the  throughput  from  1.09  to  2.22,  or  104%.  The 
absolute  increase  is  almost  constant  with  respect  to  N ,  so  for  larger  N,  the  percentage 
increase  becomes  small.  Nevertheless,  even  though  the  increase  in  maximum  throughput 
may  be  small,  by  using  channel  load  sensing,  the  system  can  be  made  more  stable.  We  note 
that  in  Fig.  6.5,  the  throughput  varies  only  slightly  over  a  wide  range  of  offered  traffic. 

For  the  model  of  non-zero  propagation  delay,  we  plot  the  normalized  throughput  S/N 
vs.  l/h  for  the  half-duplex  model  with  Af  =  12.  Results  are  shown  in  Fig.  6.7  for  N  =  64 
and  for  N  =  128  chips  per  bit,  which  give  cut-offs  of  L  =  7  and  L  —  13  respectively. 
These  results  were  calculated  for  no  channel  load  sensing  and  for  K  at  the  optimum 
channel  load  sense  point.  It  is  seen  that  even  for  no  channel  load  sensing,  the  throughput 
decreases  as  l/h  increases,  due  to  the  users  being  busy  during  the  holding  periods.  Also, 
as  expected,  the  significant  improvement  from  channel  load  sensing  for  N  =  64  degrades 
as  l/h  increases,  and  there  is  almost  no  improvement  for  l/h  =  1.0. 

7.  Conclusions 

In  this  report,  we  developed  the  Markovian  model  of  a  packet  radio  network  for  several 
specific  fully  connected  topologies.  We  Produced  a  refinement  which  allows  general  bit 
error  rate  functions  and  incorporated  a  model  of  a  DS-BPSK  CDMA  radio  channel  with 
convolutional  FEC  coding.  Throughput  and  normalized  throughput  were  found  for  various 
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Normalized  Throughput 


Figure  6.5.  Throughput  5  versus  offered  traffic  for  several  values  of  K 
for  an  infinite  population  with  channel  load  sensing. 


Throughput 


Figure  6.6.  Throughput  5  vs.  channel  load  sense  point  K 
for  an  infinite  population. 
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Normalized  Throughput 


Figure  6.7.  Normalized  throughput  S/N  versus  delay  parameter  l/h 
for  the  half-duplex  model  with  channel  load  sensing. 
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values  of  Eb/No  and  chips  per  bit  N.  Next,  we  accounted  for  the  effect  of  both  the  receivers 
and  the  transmitters  belonging  to  the  finite  user  population,  as  would  be  the  case  for  a 
half-duplex  radio  network.  This  had  a  large  impact  on  throughput.  We  then  considered 
a  channel  load  sensing  protocol,  and  found  that  this  resulted  in  increased  throughput. 
Finally,  a  model  of  non-zero  propagation  delay  was  developed.  We  showed  that  increasing 
the  delay  degrades  the  improvement  due  to  channel  load  sensing,  a  result  familiar  from 
CSMA  analyses. 


Appendix  A.  Throughput  Equation  For  Finite  Transmitters 

To  find  the  throughput  for  a  user  t,  we  expand  the  state  space  X(t)  to  explicitly 
indicate  the  state  of  user  »,  giving  a  new  state  space  £*(t) 

We  define  the  following: 

j*  is  the  number  of  users  transmitting  not  including  user  i 

c,  is  the  state  of  user  t,  with  c,  E  {idle,  busy},  where  t  is  busy  when  transmitting 

5*  is  the  expanded  state  space,  5*  =  {(;*,  c,)  :  0  <  j*  <  M  —  1,  c,  E  {idle,  busy}} 

jr|y»  c  j  is  the  steady  state  probability  of  being  in  state  (j* ,  c, ) 

is  the  fraction  of  time  that  user  i  is  successfully  transmitting  packets  which 
found  the  network  in  state  (/*,  idle)  upon  arrival. 

The  resulting  Markov  chain  is  shown  in  Fig.  Al. 

The  analysis  then  follows  exactly  the  solution  to  a  more  general  case  for  multihop 
networks  found  by  Brazio  and  Tobagi  [BRAZ84].  First,  we  note  that  the  times  of  successive 
transitions  from  the  state  (/*,  idle)  to  the  state  (;*,  busy)  are  regeneration  points  for 
the  Markov  chain  Consider  the  cycles  defined  by  the  time  intervals  between  two 

successive  regeneration  points.  Let  Ck(j*,i)  denote  the  length  of  the  fcth  cycle,  and  let 
Tk(j*,i)  denote  the  successful  length  of  the  packet  whose  arrival  initiated  cycle  k.  The 
sequence 

is  a  sequence  of  i.i.d.  pairs  of  random  variables.  Let  E(Ck{j* ,t))  and  E(Tk{j* ,t))  denote 
the  average  values  of  Ck{j*,i )  and  Tk{j*,i)  respectively.  From  the  theory  of  renewal 
processes, 

£(;'*  >*)  =  •’  l)j  with  Probability  1  (A1) 
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(M  -  1)A  (M  -  2)A 


Figure  Al.  Expanded  state-transition-rate  diagram 


for  the  finite  transmitters  model. 
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Also,  as  shown  by  Brazio  and  Tobagi  in  [BRAZ84J, 

E(CtU\>))  =  ™ -  (A2) 

*  l/,  idle) 

If  user  i  is  idle,  and  there  are  j*  other  radios  transmitting,  the  total  number  of  users 
transmitting  is  j*.  Thus,  the  event  {^*  =  ( j  * ,  idle)}  is  the  same  as  the  event  {£  = 
j  and  t  is  idle)},  when  j*  =  j.  In  other  words, 


idle)  vjP*iDLEU)  j  ~  J  (A3) 

so  Eqn.  A2  becomes 


E{Ck{j\i )) 


1 

^  n]PtiDLi;(i) 


(A4) 


Now,  because  all  users  are  identical,  the  expected  successful  length  of  a  packet  that  is 
transmitted  given  j*  other  radios  transmitting  when  it  arrived  does  not  depend  on  which 
particular  radio  i  the  packet  arrived  at.  Consequently,  E[Tk{j*,i))  is  simply  E(T„\j)t  for 
J*  =  j- 

=  *jP,IDLEU)E(Tt\j))  j*  =  j  (A5) 

The  total  throughput  is  simply  the  sum  of  conditional  throughputs,  so 

Si=  eW.O  (A6) 

>•=0 

=  Z  *jPi,DLEV)M{T.lj)  (A7) 

7=0 


Finally,  we  notice  that  if  L  <  M,  S(j*,i)  =  0  for  j*  >  L,  so  less  computation  is  required 
for  smaller  L. 
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Appendix  B.  Solution  to  E(Tt\j) 


This  derivation  of  E(T,\j)  from  the  matrix  B.  is  also  from  [BRAZ84]. 

P*(r)  is  the  probability  transition  matrix  of  the  auxiliary  Markov  chain,  defined  as 

p;-y ,(t)  =  Pr(r(r)  =  j'ir(o)  =  j),  j  e  $««x,  j'  e  Saux 


Similar  to  the  corresponding  rate  transition  matrix  R*,  P*(r)  has  the  form 


P(r)  Ps(r)  P  F(r) 
P*(r)  =  I  0  1  0 

0  0  1 


Because  of  this  structure,  the  forward  Kolmogorov  equations 

d 


dr 


P‘(r)=P-(r).R*  P*(0)  =  I 


can  be  broken  down  to 


d 

dr 

d 

dr 

d 

dr 


P(r)=P(r)R 
Ps(r)  =  /iP(r)l 
PF(r)  =  P  (r)<p 


P(0)  =  I 

p5(0)=0 
PF(0)  =  0 


which  has  the  solution 


(Bl) 


(B2) 


(B3) 


P(t)  -  eRr 

r  >  0 

Ps(r)  =  ^(eRT-I)R-1l 

T  >  0 

(B4) 

Pf(0  =  (^Rr-i)B”V 

T  >  0 

Since  all  the  states  except  Success  and  Failure  are  transient,  eRr  -»  0  as  r  oo.  By 
construction,  the  j  +  1st  element  of  the  column  vector  P^Coo)  is  the  probability  that  a 
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transmitted  packet  which  found  the  network  in  state  j  on  arrival  will  be  absorbed  in  the 
state  Success,  which  we  defined  in  Section  3. A  as  Ps\j-  Thus,  from  Eqn.  B4,  we  find 
Ps(oo)  =  fi(— I)E_11  =  — /xE-1l.  Ps(oo)  is  equal  to  9  from  Section  3. A. 

Now,  let  T  be  a  random  vector  giving  the  time  to  absorption  in  Success.  We  construct 
T  as  column  vector  with  rows  indexed  by  the  non-absorbing  states  in  $avz  in  the  same 
order  as  the  rows  of  E.  The  j  4-  1st  element,  [T]y+j ,  is  the  random  variable  which  is 
the  'successful  length  of  a  packet  which  arrived  to  find  the  network  in  state  j  and  was 
transmitted. 

£  Pr(|T|,+i)  =  PsM  r>0  (BS) 

and 

E(T|T  <  oo)  =  jf “  (Ps(oo)  -  Ps(r))  dr 

=  -ft  I™  eRTR~ll  dr  =  hR~21  (B6) 

J  0 

The  desired  quantity  E(T,[j)  is  simply  the  j  +  1st  element  of  this  vector,  so 

E(T, |j)  =  J5([T|J+,  |  [T]J+I  <  oo)  =  [^R-21]J+1  (B7) 
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Appendix  C.  Viterbi  Decoder  performance 

In  this  appendix,  we  summarize  the  analysis  required  to  derive  the  performance  of  the 
decoder.  The  analysis  is  due  mainly  to  Viterbi  [VITE71],  Also,  an  excellent  treatment 
of  convolutional  coding  and  of  the  Viterbi  decoding  algorithm  is  presented  by  Clark  and 
Cain  in  [CLAR81]. 

C.l  Characteristics  of  Convolutional  Codes 

Unlike  block  codes,  convolutional  codes  do  not  require  a  frame  structure.  Instead,  the 
output  symbols  from  the  encoder  at  any  time  are  a  linear  combination  (modulo  2)  of  the 
previous  k  bits,  so  the  output  stream  depends  on  the  input  bits  in  a  sliding  window  fashion. 
For  a  rate  m/n  code,  n  symbols  are  output  for  each  m  bits.  Typically,  for  m  =  1  codes,  a 
shift  register  will  store  the  last  k  bits,  and  each  of  the  n  symbols  will  be  generated  by  a 
different  linear  binary  function  of  these  bits. 

The  major  parameters  of  the  code  are  its  constraint  length  k  and  its  rate  m/n.  We 
only  consider  a  commonly  used  rate  1/2  constraint  length  7  code  in  the  network  model, 
but  for  illustrative  purposes  we  will  also  discuss  a  rate  1/2  constraint  length  3  code. 

Convolutional  codes  are  group  codes,  meaning  that  a  set  of  symbols  of  a  given  length 
form  a  group.  One  convenient  property  of  group  codes  is  that  for  a  binary  symmetric 
channel  (BSC),  *he  probability  of  error  does  not  depend  on  the  actual  data  stream  being 
sent.  Thus,  we  can  analyze  the  probability  of  error  for  the  data  stream  of  all  zeros  and 
this  will  yield  the  probability  of  error  for  any  data  stream.  Furthermore,  because  the 
codes  are  linear,  the  all  zero  data  stream  will  result  in  all  output  symbols  being  zero.  For 
convenience,  we  denote  the  symbol  stream  of  m£  zeros  by  C(1  .  Also,  without  loss  of 
generality,  we  assume  that  the  initial  loading  of  the  encoder  shift  register  is  a  sequence  of 
k  -  1  zeros. 
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C.2  Performance  Analysis  [VITE71] 


Two  performance  measures  are  commonly  derived  for  evaluating  codes.  The  most 
widely  used  measure  is  Pg,  the  probability  of  bit  error.  This  is  the  long  term  proportion 
of  bit  errors  in  a  very  long  sequence  of  bits.  Another  measure  is  Pe ,  the  probability  of 
first  error.  This  is  the  probability  of  a  bit  error  at  the  output  of  the  decoder,  given  that 
the  first  A:  —  1  bits  were  known  by  the  decoder,  and  given  that  no  errors  have  occurred 
up  to  the  present.  In  the  network  model,  we  assume  that  a  known  header  precedes  the 
packet  transmission.  This  header  will  allow  the  decoder  to  be  initialized  to  the  correct 
state.  Ouce  initialized,  we  are  interested  in  the  first  occurrence  of  an  error  in  the  decoder 
output  bit  stream,  as  this  will  cause  the  entire  packet  to  be  declared  in  error. 

The  Viterbi  decoding  algorithm  is  a  maximum  likelihood  decision  rule,  in  which  the 
codeword  which  is  closest  in  Hamming  distance  to  the  received  codeword  is  chosen  as  the 
estimate  of  the  transmitted  symbol  stream,  and  the  corresponding  bit  stream  is  output 
as  the  estimate  of  the  source’s  information  stream.  We  denote  the  Hamming  distance 
between  Cn  and  Cj  by  du[Cn,  C  ?).  The  probability  of  error  given  that  Cq^  was  sent  is  the 
probability  that  the  first  m£  received  symbols,  is  closer  to  some  other  codeword  of 

/  p  \ 

length  £  than  it  is  to  C0  .  This  is  the  union  of  events, 

<W*l£l,Ci£l)} 

0=1 

/  p  \  i  n  \ 

In  the  case  of  a  tie  between  Cn  and  Cq  ,  one  is  chosen  at  random.  We  can  bound  the 
probability  of  this  event  with  the  union  bound,  so  that 

2£-l 

Pr(error)  <  £  Pr(d„(Rlc> ,Cln)  <  dn(K,c> 

0=1 

+  1/2  £Pr(dH(£l£l,C!,£))  =  dgfX^.C1,"))  (Cl) 

0  =  1 
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This  union  bound  is  the  most  common  method  used  for  analyzing  decoder  performance. 

We  can  model  the  encoder  as  a  finite  state  machine.  Transitions  between  states  occur 
at  each  new  input  bit.  Because  the  oldest  bit  of  the  k  bits  in  the  shift  register  does  not 
affect  the  next  state,  we  only  need  2k~l  states.  The  finite  state  machine  representation  of 
the  rate  1/2  k  —  3  encoder  is  shown  in  Fig.  Cl. 

There  is  a  one-to-one  mapping  between  codewords  and  paths  through  the  finite  state 

(L\ 

machine.  Consider  two  paths  and  of  length  £,  corresponding  to  codewords  C[a 
and  C[f\  which  both  start  in  the  state  zero  and  end  in  the  same  state.  Let  be  the 
set  of  all  codewords  of  length  Ci  >  L  for  which  the  first  t  bits  are  the  codeword  Ca  . 
Every  codeword  C^1*  G  A is  the  concatenation  of  C'a'1  and  some  sequence  C^2\  where 
C  +  £2  =  C\.  Similarly,  we  form  the  set  B from  C^:f  \  and  find  that  G  B is 
the  concatenation  of  C*/*  and  Also,  we  denote  the  last  mCi  symbols  of  by 

We  wish  to  find  the  sequences  in  and  at  minimum  Hamming  distance  from 
C^1.  We  note  that 


dH(^c'\C\c'])^dn{Z{C\ClP)  +  dn(Z^\c[Ci)) 


dH(Zlc']XiC'])  -dn{^\C(P)  +  dn{R^\C[^])  (C3) 

Let  G  be  the  codeword  which  minimizes  djj{Z^'  *,  C^-f'  *),  and  similarly  for 

C^''>  G  .  Because  the  starting  state  is  the  same,  C [52*  and  range  over  the  same 

sequences.  This  implies  that  dg(Z[^2 must  equal  dg{Z^I2\ Therefore, 
the  codeword  which  is  uniquely  closest  to  £1^1’  will  be  in  A^‘  1  only  if  dg(Z{C\ C^1)  < 
and  will  be  in  B only  if  dtf(Z^^ ,  C!f^ )  >  dg(Z^\C^^).  In  the  case 
of  a  tie,  dj{{Z^) ,  C\P )  =  d{j(Z^\C{f]),  the  codewords  C^''1  and  will  have  equal 


l£\ 

Hamming  distances  to  C0  ,  so  one  must  be  chosen  at  random.  We  can  make  this  random 
choice  at  bit  £  without  increasing  the  probability  of  error.  Therefore,  no  matter  how  the 
paths  compare,  at  bit  £  we  can  eliminate  one  of  the  paths  entering  a  given  state.  In  fact, 
we  can  eliminate  all  but  one  path  entering  every  state.  We  call  the  remaining  paths  the 
survivors. 


Because  of  the  known  starting  state,  during  the  first  k  -  2  bits,  there  are  a  limited 
number  of  states  that  the  encoder  can  be  in.  For  example,  it  cannot  be  in  the  all  ones 
state  since  all  of  the  zeros  will  not  have  been  shifted  out.  After  A:  —  1  bits,  it  can  be  in  any 
state,  so  the  bit  can  be  mapped  into  any  of  the  2fc  transitions.  There  are  at  most  2k~l 
survivors.  At  each  bit  a  comparison  must  be  made  between  the  two  paths  entering  every 
state,  or  at  most  2,;  comparisons.  The  decoder  must  retain  the  2fc_1  survivors  and  the 
Hamming  distance  between  the  received  symbol  sequence  and  each  of  the  2k~l  survivors. 

For  a  finite  length  message,  a  known  sequence  of  k  -  1  bits  must  be  appended.  This 
will  give  a  uniquely  specified  ending  state,  so  the  decoder  will  be  able  to  decide  between 
the  2/:  “ 1  survivors. 


( £. ) 

Assuming  that  the*  all  zero  codeword  Cu  was  sent,  an  error  event  occurs  at  bit  £  when 
the  path  merging  with  the  state  zero  is  chosen  as  a  survivor,  eliminating  C0  .  For  a  fixed 
length  message,  if  we  also  assume  that  the  known  concluding  bits  are  a  sequence  of  k  —  1 
zeros,  the  decoder  will  always  ch<  >ose  a  codeword  which  ends  in  the  state  zero,  so  again 
the  only  ."mnparisons  we  need  to  consider  are  those  between  Cl  and  the  path  entering 


the  state  zero.  Denote  the  codeword  corresponding  to  the  path  entering  the  state  zero  at 
bit  £  by  Cl£)  ■  The  difference  in  Hamming  distance  between  and  and  ^ '  and 
C;1"1  only  depends  or.  th  ?e  symbols  of  Z11-1  for  which  C ,J’1  -  1.  Let  i  be  the  number  of 


symbols  for  which  C {  1.  The  probability  that  C'^'  is  chosen  as  the  survivor  is  P,,  the 
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probability  that  more  than  half  of  these  t  symbols  are  in  error,  which  is 


p,  =  { 


£ 

e=(*  +  l)/2 

>(•' 

2  Vi/2, 


P'(  1  -  P)-’, 


i  odd; 


(C4) 


P*/2(i  -Py/2+  £  * 


even. 


e=t/2+ 1 


where  p  =  Pe , yt„M,  the  channel  symbol  error  rate. 

We  wish  to  find  the  number  of  codewords  which  begin  in  state  zero,  return  to  state  zero 
for  the  first  time  at  bit  £  and  have  i  ones.  This  is  done  by  modifying  the  state  diagram  as 
follows.  Split  the  zero  node  into  a  source  node  and  a  sink  node.  Assign  a  gain  of  LN^D1 
to  the  transitions,  where  /?  is  the  weight  of  the  input  bit  (i.e.,  0  for  a  0,  1  for  a  1),  and  7 
is  the  weight  of  the  output  bit  (i.e.,  0  for  00,  1  for  01  or  10,  and  2  for  11).  For  a  path  of 
length  a1  with  ft1  ones  in  the  input  bits  and  7'  ones  in  the  output  symbols,  the  product  of 
the  transition  gains  is  La' N  :1' D^' .  The  modified  state  diagram  is  given  in  Fig.  C2. 

If  we  sum  together  all  possible  paths  from  the  source  to  the  sink,  we  will  find  the 
overall  gain  of  the  network,  T(D,  L,  TV).  For  example,  for  the  rate  1/2  k  =  3  code,  this  is 


T{D,  L,  TV )  =  DbLzN  +  D*LA{  1  +  L)N 2  +  ...  +  Db+}Lz+}{  1  +  L)J  Nl+]  +  ..  (C5) 

This  gain  can  be  found  directly  from  the  modified  state  diagram  by  solving  the  set  of  linear 
equations  giving  the  dependencies  between  the  nodes.  For  example,  for  the  rate  1/2  k  =  3 
code,  we  find  T(D ,  L,  N )  directly  as 


T(D,L,N ) 


D°LZN 

1  -  DL{  1  +  L)N 


(C6) 


Since  we  are  not  interested  in  the  number  of  ones  in  the  input  bits,  set  TV  =  1  to  get 


T(D,L) 


DbLz 

1  -  DL{  1  +  L) 


(C7) 
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Figure  C2.  Modified  state  diagram. 
Rate  1/2  k  =  3  code. 
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We  define  a,  to  be  the  sum  of  the  coefficients  of  all  terms  of  T(D,  L)  for  which  the 
exponent  of  D  is  i  and  the  exponent  of  L  is  less  than  or  equal  to  C.  is  the  number 
of  codewords  of  length  less  than  or  equal  to  C,  with  t  ones.  The  union  bound  on  the 
probability  that  an  error  event  occurs  at  bit  C  is 

m  £ 

(C8) 

«=1 

This  is  overbounded  by  letting  C  go  to  infinity.  We  define  o,  =  lim£_00a-^>  so 

OO 

Pe  <  £  a>Pi  (C9) 

«=1 

for  every  bit.  We  find  the  coefficients  a,  by  evaluating  T(D,L)  at  L  =  1,  since 

T(D,L)\L=l  =  f2atD'  (CIO) 

»  =  i 

An  additional  simplification  is  realized  by  using  the  bound  Pt  <  (2^/p(l  —  p))‘.  This 
yields 

Pe  <  t  *Vfyl-r)Y  =  T(D) (01 1) 

«  =  1 

However,  this  is  a  very  loose  bound  for  large  values  of  p.  Because  we  are  interested  in  the 
decoder  performance  up  to  about  10— 2 ,  it  is  necessary  to  use  the  more  exact  expression  of 
equation  C9. 

C.3  Numerical  evaluation  of  Decoder  Performance 

Several  steps  were  required  in  order  to  derive  numerical  results  for  the  performance. 
First,  the  transfer  function  T(D)  had  to  be  calculated.  Next,  the  expression  had  to  be 
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stated  in  the  form  of  a  power  series  expansion  to  find  the  coefficients  a,.  Finally,  the  terms 
Pi  had  to  be  calculated  and  the  summation  of  Eqn.  C9  evaluated. 

A  computer  program  was  written  to  give  the  equations  relating  the  states  of  the  mod¬ 
ified  state  diagram  for  the  rate  1/2  constraint  length  7  code  being  modeled.  The  MAC- 
SYMA  (©1976,1983  Massachusetts  Institute  of  Technology,  ©1983  Symbolics,  Inc.)  sym¬ 
bolic  manipulation  program  was  then  used  to  solve  the  64  simultaneous  equations,  giving 
an  expression  for  T(D)  as  the  ratio  of  two  polynomials.  This  expression  is  given  in  Table 
Cl. 

Table  Cl.  Polynomial  Ratio  Solution  to  T(D) 

Numerator  -  D80  +  16D76  -  120D72  -  D70  +  562D68  +  8 D66  - 1838 £>64  -  20 D62  +  44  29 D60  - 
18£>58  -  80  68D56  +235£>54  + 1 12 18D52  -  678  Z?50  —  1 1900Z>48+  1097Z?46+9575Z)44-  1094Z?42- 
5841D40  +  61  ID38  +  2795 £>36  -  49D34  -  1156D32  -  22  8  D30  +  417D28  +  243D26  -  76£>24  - 
176 D22  -  15 Dw  +  93D18  +  Dlt  -  25 Du  -  6 Dn  +  llD10 

Denominator=2D  76  -  32D72  +  240  D68  +  3D66  -  1123Z?64  -  30D62  +  36  62 Z>60  +  131Z>58  - 
8766D56  -  33 ID54  +  15  7  63  D52  +  561D50  -  21403D48  -  782 D46  +  21746D44  +  1184  D42  - 
16133Z?40  -  1960Z?38  +  83  44  D36  +  2  8  07 D34  -  2751  £>32  -  3064Z?30  +  389I>28  +  2509D26  + 
267 D24  -  1601 D22  -  403  D20  4-  844D18  +  262D16  -  345D14  -  81 D12  +  85 D10  +  40 D8  -  30D6  - 
6  D*  -  4  D2  +  1 

The  denominator  was  then  expanded  using  the  identity 

—  1  +  x  -f-  X2  -f-  x3  +  . . .  (C12) 

1  -  x 

The  expansion  was  calculated  out  to  64  terms.  This  was  then  multiplied  by  the  numerator, 
and  the  first  64  terms  were  retained.  The  coefficients  are  listed  in  Table  C2. 

We  note  that  for  large  i,  the  number  of  codewords  of  Hamming  weight  i  given  by  a,  is 
very  large.  For  a  low  probability  of  symbol  error  P>  will  be  small  enough  to  offset 
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Table  C2.  Coefficients  o. 


1 

a. 

i 

ai 

t 

10 

11 

54 

2.912797332  x 

1017 

98 

1.243940152  x  1034 

12 

38 

56 

1.660510362  x 

1018 

100 

7.091380930  x  1034 

14 

193 

58 

9.466139591  x 

1018 

102 

4.042612776  x  1035 

16 

1331 

60 

5.396401324  x 

1019 

104 

2.304588936  x  1036 

18 

7275 

62 

3.076348764  x 

102° 

106 

1.313786523  x  1037 

20 

40406 

64 

1.753746820  x 

1021 

108 

7.489557000  x  1037 

22 

234969 

66 

9.997656840  x 

1021 

110 

4.269602640  x  1038 

24 

1337714 

68 

5.699405471  x 

1022 

112 

2.433989993  x  1039 

26 

7594819 

70 

3.249083581  x 

1023 

114 

1.387554717  x  104° 

28 

43375588 

72 

1.852218477  x 

1024 

116 

7.910090375  x  104° 

30 

247339453 

74 

1.055901826  x 

1025 

118 

4.509337848  x  1041 

32 

1409277901 

76 

6.019423091  x 

1025 

120 

2.570656827  x  1042 

34 

8034996288 

78 

3.431517347  x 

1026 

122 

1.465464930  x  1043 

36 

4.580875611  x  1010 

80 

1.956219246  x 

1027 

124 

8.354236400  x  1043 

38 

2.611287754  x  10u 

82 

1.115189974  x 

1028 

126 

4.762534003  x  1044 

40 

1.488634502  x  1012 

84 

6.357409478  x 

1028 

128 

2.714997404  x  1045 

42 

8.486419243  x  1012 

86 

3.624194641  x 

1029 

130 

1.547749762  x  1046 

44 

4.837861791  x  1013 

88 

2.066059586  x 

103° 

132 

8.823320927  x  1046 

46 

2.757937903  x  1014 

90 

1.177807109  x 

1031 

134 

5.029946958  x  1047 

48 

1  572231420  x  1015 

92 

6.714373564  x 

1031 

136 

2.867442606  x  1048 

50 

8.962880896  x  1015 

94 

3.827690629  x 

1032 

52 

5.109505443  x  1016 

96 

2.182067383  x 

1033 
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the  coefficient  a,-,  so  only  the  first  few  terms  of  the  summation  are  significant.  However, 
for  Pe^ymbcl  =  0.04715,  which  corresponds  to  PE  =  10-2,  the  term  al3SPm  is  equal  to 
2.16  x  10-4,  or  about  2%  of  the  total.  For  Pe,»ymbol  >  0.04715,  the  power  series  yields 
a  PE  >  10~2.  However,  in  this  range,  terms  of  the  summation  of  order  higher  than  138 
become  significant.  Because  the  numerical  calculation  truncates  these  terms,  the  power 
series  result  is  not  the  true  union  bound.  Of  course,  the  bound  becomes  quite  loose,  so 
the  partial  sum  is  probably  greater  than  the  actual  PE.  Nevertheless,  we  use  the  looser 
bound  PE  <  1  for  Pe,Kymbol  >  0.04715. 
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Appendix  D.  Probability  of  Correct  Packet  Reception 

D.l  Fixed  Length  Packets  [TAEP84] 

For  a  packet  of  length  L  being  sent  over  a  binary  symmetric  channel,  the  probability  of 
the  packet  being  correctly  decoded,  Pc ,  was  bounded  by  Taiple  [TAIP84]  as  follows.  For 
correct  decoding,  we  assume  that  the  starting  state  is  known  and  that  a  known  tailer  of 
k  —  1  bits  is  transmitted,  so  a  total  of  i  -f  k  -  1  bits  are  encoded,  and  m(£  -f  k  —  1)  symbols 
are  transmitted.  Again,  we  choose  the  starting  state  and  the  tailer  to  be  A:  —  1  zeros.  There 
are  a  total  of  2£  codewords  which  are  candidates  for  comparison  to  the  received  sequence, 
since  the  decoder  will  only  choose  a  codeword  which  starts  and  stops  in  the  all  zero  ^iate. 
We  index  the  paths  corresponding  to  these  codewords  as  £,.  Let  the  event  E,  be  the  event 
that  path  £,  is  not  chosen  by  the  decoder.  The  packet  will  be  correct  if  the  path  0  is  chosen 
by  the  decoder,  or  equivalently,  if  E,  does  not  occur  for  any  i  >  0.  Thus, 

2£-l 

Pc  =  Pr{  n  E.)  (Dl) 

t=l 

Taiple  [TAIP84]  proves  that  for  the  events  E,, 

2£-l 

Pc  >  n  Pr(E.)  (D2) 

t=i 

We  can  group  the  non-zero  codewords  according  to  the  first  time  at  which  they  return 
to  the  zero  state.  Let  Bx  be  the  set  {t|£,  returns  to  the  zero  state  at  bit  x}.  The  soonest 
a  non-zero  codeword  can  return  to  the  zero  state  is  in  k  bits,  since  a  one  followed  by  all 
zeros  still  must  pass  through  all  k  stages  of  the  encoder  shift  register.  Thus,  the  first 
non-empty  set  Bx  is  Bk .  The  sets  {Bx\k  <  x  <  C  +  k  —  1}  form  a  partition  of  the  indices 
{1,2,...,2£-1},so 

2£-l  C+k-l 

Pc>  n  pr(E.)  =  n  n  jme.)  (ao 

i  =  I  x  =  t :  ieDi 
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Let  {£,'}  be  the  set  of  paths  starting  at  bit  —  oo,  passing  through  any  arbitrary  state 
at  bit  zero  and  ending  in  the  zero  state  at  bit  C  +  k  —  1.  We  can  include  in  the  product 
Pr(E,)  the  terms  corresponding  to  paths  which  first  return  to  the  zero  state  at  bit 
x  but  which  did  not  start  in  the  initial  state  of  all  zeros  at  bit  0  ,  and  still  retain  the 
inequality  of  D3,  since  the  probabilities  of  the  extra  terms  are  all  less  than  1  Let  B'z 
be  the  set  {i|£'  returns  to  the  zero  state  at  bit  x}.  The  set  of  all  paths  {£'|i  e  B'x}  are 
equivalent  for  all  x,  so  we  let  B1  =  B'x.  The  inequality  of  D3  now  becomes 


Pc  >  (  II  Pr(E.) 

yieB' 

f 

>  I  n  (1  -  Pr(E,)) 

ieB' 


(D4) 


where  E,  is  the  compliment  of  E,  We  multiply  out  the  term  Il,e2?'(l  —  Pr(E,))  to  get  an 
infinite  sum,  which  is  lower  bounded  by  the  sum  of  lower  order  terms  1  -  EieB'  Pr(Et). 
However,  the  event  E,  is  the  event  that  path  is  chosen  as  the  correct  path,  and  this  is 
summed  over  all  paths  which  are  returning  to  the  zero  state.  This  sum  is  exactly  the  sum 
found  in  the  union  bound  earlier.  Thus, 

Pc>  (l-  E  pr(E,))  =(l-/’£)<:  (D5) 

V  ieB'  J 


D.2  Varying  Symbol  Error  Rate 

Using  a  similar  approach,  we  can  find  an  approximation  to  a  bound  for  the  probability 
of  correct  packet  reception  for  a  variable  length  packet  being  transmitted  on  a  channel 
with  varying  symbol  error  rate,  as  in  the  network  model.  This  approximation  results  in 
a  model  of  decoder  performance  which  is  memoryless,  and  therefore  can  be  incorporated 
into  the  Markovian  network  model. 
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The  analysis  of  the  Viterbi  decoding  algorithm  assumes  that  infinite  decoder  memory 
is  available,  so  that  the  chosen  codeword  will  be  the  one  which  is  closest  to  the  received 
symbol  stream.  In  a  practical  implementation,  a  finite  memory  must  be  chosen,  and  a 
sub-optimum  choice  must  be  made  when  the  memory  overflows.  However,  the  memory 
can  be  made  large  enough  that  the  probability  of  this  occurring  is  negligible.  Simulations 
have  shown  that  a  memory  of  approximately  5 (k  -  1)  information  bits  is  sufficient  for  a 
rate  1/2,  constraint  length  k  decoder  [CLAR81].  In  the  case  of  the  k  —  7  code,  this  is  30 
bits. 

Even  if  the  packet  length  is  not  fixed,  for  a  given  packet  of  length  £  bits,  the  inequality 
of  D3  is  still  valid.  We  again  include  those  paths  corresponding  to  codewords  of  length 
greater  than  x  in  the  product,  to  give 

£  +  k- 1  £  +  k-l 

pc  >  n  n  pr(E.)  =  n  n  a  -  pr^n  (w) 

r-k  ieB'j  x~k  I  eB'j 

Also,  the  inner  product  can  be  lower  bounded  by  the  sum  of  lower  order  terms,  so 

11(1-  Pr(E,))  >  1  -  E  Pr(E.)  (D7) 

i eD'j  ieD’x 

We  now  partition  B x  into  +  B'Z,  those  paths  which  diverged  from  the  all  zeros  state 
within  the  last  30  bits,  and  ~  B'x,  those  paths  which  diverged  more  than  30  bits  ago.  The 
decode  memory  is  chosen  such  that 

E  Pr(E.)«  E  Pr(E.)  (D8) 

ie-D'x  i €+B'j 

Thus,  we  approximate  the  sum  £lGz?'  Pr(E.)  by  E,e+z?i  Pr(E«).  To  find  Pr(E,)  for  i  e  +B'X, 
we  must  know  the  probability  of  symbol  error  during  each  of  the  previous  30  bits.  However, 
the  sum  is  monotonic  in  Pe,*ymM-  We  can  bound  the  sum  by  taking  the  nighest  and  lowest 
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values  of  Pe, symbol  occurring  in  the  last  30  bits  and  finding  Pr(E.)  if  all  symbols  had  the 
extreme  value  of  Pe, symbol- 

We  approximate  the  performance  by  using  the  current  value  of  Pe^ymbol  as  the  value 
experienced  by  the  last  60  symbols.  By  u&Ing  the  bounds  on  Pe,lymbol,  we  show  in  Section 
4.C  that  this  approximation  is  valid. 

If  the  symbol  error  rate  for  all  60  symbols  is  Pe,symbol[i)i  we  can  bound  the  sum  by 

E  Pr(Bi)  <  P£(j)  (D9) 

ie+  J', 

where  Pe(j  )  is  the  probability  of  first  error  given  a  symbol  error  rate  of  Pe,,ymboi{j )•  If  U 
is  the  state  at  bit  x,  using  the  memoryless  approximation  we  find 

C+k~  1 

Pc>  n  (l-JW*))  (D1°) 

i=k 

Let  Cj  be  the  number  of  bits  for  which  the  state  is  j. 

PciSll-PlIi))'1  (DU) 

}=  1 

As  described  in  Appendix  C,  the  decoder  model  overbounds  Pe  by  Pe  <  1  for  those 
values  of  Pe  greater  than  10~2.  We  use  L  to  denote  the  state  j  such  that  Pe{L)  <  10-2 
and  Pe{L  +  1)  >  10-2.  We  bound  P<;  by  zero  if  the  state  L  +  1  is  reached.  Thus,  we  only 
need  to  find  P( ;  for  those  cases  when  i}  —  0  for  j  >  L.  We  then  find 

Pc>  n  (D12) 

j=i 
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Appendix  E.  Throughput  Equation  for  the  Half-duplex  Model 


As  in  Appendix  A,  we  expand  the  state  space  to  explicitly  indicate  the  state  of  user  i 
We  define  the  following: 

t*  is  the  number  of  users  transmitting  not  including  user  t 
r*  is  the  number  of  users  receiving  not  including  user  i 


c,  is  the  state  of  user  », 


S *  is  the  expanded  state  space, 


i  is  idle 

i  is  transmitting 
i  is  receiving 


0  <  f*  <  M  -  1  ) 


S' 


(t',r‘,c,)  : 


0  <  r*  <  t*  +  c, 
r  +  r*  <  M  -  1 


c,  €  {-1,0,1} 


r.  p  j  is  the  steady  state  probability  of  being  in  state  (t*,r*,c,) 

5(<*,r*,t)  is  the  fraction  of  time  that  user  t  is  successfully  transmitting  packets  which 
found  the  network  in  state  (<*,r*,0)  upon  arrival. 

We  again  have  Cfc(<*,r*,t)  and  7\(<*,r*,t)  as  in  Appendix  A. 

£(Ct(‘V,.))=  (El) 

The  event  {X *  =  (<*,r*,0)}  is  the  same  as  the  event  {X  =  (<,r)  and  t  is  idle)},  t  =  t*  and 
r  =  r*,  so 

♦  ,,  >  t  =  t*  and  . 

7r(<*.r,.0)  —  K(t.r)Psi{t,  r)  r  -  r*  (®»2) 

80 


and 


E(Ck{t',r',x)) 


1 

**U,r)Psi\t,r) 


t  —  t*  and 
r  =  r* 


(E3) 


By  definition,  the  quantity  E(T,\(t,r))  is  conditioned  on  the  packet  capturing  the 
receiver.  If  the  receiver  is  not  captured,  this  corresponds  to  an  immediate  transition  to 
the  state  Failure  in  the  auxiliary  Markov  chain.  Since  ,  r* ,i)  is  not  conditioned  on 

capture,  we  have 


E{Tk(t*  ,r*  ,i))  =  E{T„\(t,  r))  •  Pr(receiver  is  captured  jf,  r)  *  *  an<^  (EM) 

r  —  r* 

Pr(receiver  is  captured  \t,r)  is  equal  to  r) at.  Given  these,  we  can  find 


t)  = 


E(Tk{f,r* 

E{Ck{t',r' 


0) 

0) 


=  ^{l,r)Psi(t,r)PDI(t,r)aiE(T.\(t,r)) 


t  =  t*  and 
r  =  r* 


(E5) 


for  ihose  states  {(t* ,  r*)  :  (t1 ,  r*,  c,j  €  5‘,  for  any  c, }.  Since  ^'(T'gKt,  r))  =  0  for  t  >  L,  the 
throughput  can  be  found  by  summing  over  the  states  S',  so 


&i=  H  Air(,.r)Ps/(<,r)PDl(f,r)a,£'(r,Kf,r))  (E6) 

(f.r)eS' 
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Appendix  F.  Throughput  Equation  for  Delay  Models 


F.l  Receivers  Available  for  all  Transmissions 

The  derivation  is  very  similar  to  those  of  appendixes  A  and  E. 

We  define  the  following: 

w *  is  the  number  of  users  in  the  holding  state  not  including  user  i 
t *  is  the  number  of  users  transmitting  not  including  user  i 


(0 


c,  is  the  state  of  user  i,  c,  =  <  1 


12 


t  is  idle 
t  is  holding 
»  is  transmitting 


{0  <  w*  <  M-  1 
(«;*,  t*,  c,)  :  0  <  t*  <  M  —  1  —  to*  1 
0  <  c,  <  2 

7r*u..  r  )  is  the  steady  state  probability  of  being  in  state  ( w*,t*,Cj ) 

S(w*,t* ,  i)  is  the  fraction  of  time  that  user  i  is  successfully  transmitting  packets  which 
left  the  holding  state  when  the  network  was  in  state  1). 

E{Ct(w‘,r,i))=  1  (Fl) 

The  event  {X'  =  is  the  same  as  the  event  {X  =  (u>,<)  and  i  is  holding)}, 

w  =  w*  and  t  =  t‘  so 


(F2) 
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in  =  w‘  and 
t  =  <* 


I 


which  gives 


E{Ck{w*,t\i))  = 


tv  =  tv'  and 
t  -  f 


Also, 


E(Tk{tv',t',i))  =  atE(Ts\(w,t)) 


tv  —  tv'  and 
t  =  t ' 


Given  these,  we  can  find 


c,  .  E(Tk(w-,r,i )) 


=  ‘/*(w.t)pH{w>t)°‘iE(T,\(w,t)) 


tv  =  tv*  and 
t  =  V 


and 


S,  =  ]T  un[wJ)PH(tv,t)atE{T9\{tv,t)) 
(w,t)eS' 


F.2  Half-duplex  Model 

We  define  the  following: 

tv *  is  the  number  of  users  in  the  holding  state  not  including  user  t 
t*  is  the  number  of  users  transmitting  not  including  user  t 
r*  is  the  number  of  users  receiving  not  including  user  i 


(F3) 


(F4) 


(F5) 


(F6) 


c,  is  the  state  of  user  i, 


0 

1 


13 


0 


a.  =  \  1 


1-1 


t  is  idle 

i  is  holding 

t  is  transmitting 

i  is  receiving 
c,  =  0  or  1 

c,  =  2 

c,  =  3 
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0  <  to*  <  M  -  1 
0  <  t*  <  M  -  1  -  to* 

S'  =  |  (t o*,t*,r*,Cj)  :  0  <  r*  <  t*  +  a, 

to*  +  t*  +  r*  <  M  -  1 
0  <  c,  <  3 

n*w,  (,  r,  c  j  is  the  steady  state  probability  of  being  in  state  (to*,  t* ,  r* ,  c,) 

S(to*,  t* ,  r*,  i)  is  the  fraction  of  time  that  user  i  is  successfully  transmitting  packets 
which  left  the  holding  state  when  the  network  was  in  state  (to*,  t* ,  r‘,  1). 

1 


E(Ck(xv\t*,r\i))  =  — 


un: 


(F7) 


1  (w* 

The  event  { X *  =  ( xv*,t*,r *,  1)  is  the  same  as  the  event  {X  =  (w,t,r)  and  t  is  holding)} 
to  =  to*,  t  —  t *,  and  r*  =  r,  so 


*  n  (  j  \  XV  =  XV*,t  -  t*, 

and  r  =  r. 


which  gives 


E(Ck(w\t\r*,i))  = 


Also, 


£(n(to‘,<*,r*,t))  =  Pp/(to,  t,  r)or<  E{T„\{t,  r)) 


XV  =  XV*,  t  =t*, 

and  r  =  r* 


XV  =  XV* ,  t  =  **, 

and  r  —  r* 


(F8) 


(F9) 


(F10) 


Given  these,  we  can  find 

£(T*(u;VV‘,t)) 


S(to*,t*,r*,i)  = 


£(<?*(«>*, **,  »■*»*)) 

•n  —  *n*  /  —  4* 

r)PDf(wtt,  r)ott E{Tg\{xv,  t,  r))  ^  ^  ^  ’(Fll) 


and 


5,  =  £  I/,r(tr, (to,  t,  t)Pdj(xv,  t,r)aiE(Tt\(xv,  t,r))  (F12) 

(w,t,r)eS' 
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